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Introduction.  Dynamic  bifurcation  theory  in  differential  equations  is  concerned 
with  the  changes  that  occut  in  the  structure  of  the  limit  sets  of  solutions  as  para¬ 
meters  in  the  vector  field  are  varied.  For  example,  if  the  vector  field  is  the 
gradient  of  a  function  with  a  finite  number  of  critical  points,  then  the  w-limit  set 
of  each  orbit  is  an  equilibrium  point.  Thus,  one  must  be  concerned  with  how  the 
number  of  equilibrium  points  changes  with  the  parameters  (this  is  usually  called 
static  bifurcation  theory),  how  the  stability  properties  of  the  equilibrium  points 
change  and  the  manner  in  which  the  equilibrium  points  are  connected  to  each  other 
by  orbits.  If  the  vector  field  is  not  the  gradient  of  a  function,  then  other  types 
of  limiting  motions  can  occur;  for  example,  periodic  orbits,  invariant  tori,  homo- 
clinic  and  heteroclinic  orbits.  Important  questions  in  bifurcation  theory  are  con¬ 
cerned  with  the  manner  in  which  these  more  complicated  limit  sets  change  as 
parameters  vavy. 


A  person  being  introduced  for  the  first  time  to  bifurcation  theory  may  have  the 
impression  that  it  consists  of  a  collection  of  isolated  results  without  any  unifying 
principles.  Furthermore,  since  bifurcations  are  a  rare  occurrence,  perhaps  they 
could  be  avoided  if  one  were  clever  enough.  Neither  of  these  statements  are  true. 

As  an  illustration, suppose  one  ha6  a  one  parameter  family  of  vector  fields  de¬ 
pending  on  a  parameter  X,  0  <_  X  <_  1,  with  the  property  that  the  phase  portraits  of 
the  flows  for  X  *  0  and  X  =  1  are  completely  different.  Then  there  must  be  some 
point  Xq  in  (0,1)  where  the  structure  of  the  flow  changes  in  a  neighborhood  of  XQ; 
that  is,  a  bifurcation  must  occur.  This  shows  bifurcations  cannot  be  avoided.  The 
underlying  principle  in  bifurcation  theory  for  this  illustration  with  one  parameter 
familes  of  vector  fields  is  the  following.  Among  all  of  the  one  parameter  families 
Xx.  0<  X  <  1,  of  vector  fields,  characterize  those  for  which  the  bifurcations  are 
the  most  elementary.  By  most  elementary,  one  generally  means  that  a  perturbation 
of  the  one  parameter  family  will  have  the  same  type  of  bifurcations  as  the  unper¬ 
turbed  family.  This  implies  that  such  a  family  is  "transversal"  to  all  of  the 


bifurcation  surfaces  in  the  class  of  all  vector  fields.  If  the  family  X^  depends 


on  two  parameters  X 


(*1 » ^  ^  X  ^ 


1,  then  one  can  attempt  in  the  same  way  to 
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classify  those  which  are  "transversal"  to  all  the  bifurcation  surfaces.  Two  mani¬ 
folds  are  transversal  if  the  tangent  spaces  span  the  whole  space.  In  Figure  1,  we 
have  schematically  indicated  one  and  two  parameter  families  X^  which  are  trans¬ 
versal  to  the  bifurcation  surface  S. 


In  the  following,  when  we  use  the  term  codimension  one  (codimension  k)  singularity, 
we  mean  an  elementary  bifurcation  point  for  a  one  (k)  parameter  family  of  vector 
fields.  If  the  reader  keeps  this  idea  in  mind  as  he  studies  the  subject,  he  will 
recognize  that  specific  theorems  are  precise  mathematical  descriptions  of  the  above 
imprecise  remarks.  Of  course,  it  should  be  clear  that  the  same  remarks  apply  to 
vector  fields  which  depend  on  two  or  more  parameters. 

The  purpose  of  these  lectures  is  to  introduce  the  reader  to  some  of  the  basic 
ideas  in  bifurcation.  The  first  lectures  deal  with  applications  of  the  Fredholm 
alternative  and  the  method  of  Liapunov -Schmidt  to  bifurcation  near  equilibrium  and 
the  existence  of  homoclinic  orbits.  To  illustrate  the  more  global  aspects  of  the 
theory,  we  summarize  the  codimension  one  singularities  in  the  plane  and  give  some 
examples  of  codimension  two  singularities. 

k 

Throughout  the  notes  C  (X,Y)  denotes  the  set  of  functions  from  X  to  Y  which  are 

continuous  together  with  derivatives  up  through  order  k.  The  space  cJi(X.Y)  is  the 
k  b 

set  m  C  (X,Y)  with  all  derivatives  up  through  order  k  bounded  with  the  norm 

being  the  sup  of  all  derivatives  up  through  order  k. 

1.  The  Fredholm  alternative  and  Liapunov-Schmidt .  Many  problems  in  bifurcation 
theory  lead  to  the  study  of  the  zeros  of  a  function  in  the  neighborhood  of  a  given 
point.  Often,  the  analysis  consists  of  the  following  steps:  firstly,  analyze  the 
nonhomogeneous  linear  equation  (referred  to  as  the  Fredholm  alternative);  secondly, 
use  this  information  to  reduce  the  original  problem  to  one  of  lower  dimension  by 
obtaining  a  bifurcation  function  (referred  to  as  the  method  of  alternative  problems 
or  the  method  of  Liapunov-Schmidt);  thirdly,  analyze  the  bifurcation  equation; 
fourthly,  relate  the  analysis  to  dynamical  behavior. 

The  purpose  of  this  section  is  to  give  an  abstract  version  of  the  first  two 
steps. 

If  P  is  a  continuous  projection  on  any  Banach  space  X,  wc  let  Xp  denote  the 
range  of  P.  If  X,I  are  Banach  spaces,  we  let  denote  the  space  or  hounded 
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linear  operators  from  X  to  Z.  If  A  ei^X.Z),  we  let  _/f'(A)  =  {x  :  Ax  *  0},  i?( A)  = 

{z  €  Z  :  3x  6  X  3  Ax  =  z).  We  shall  also  use  the  notation  _/H[A)  =  Xy,  dl{ A)  =  Z£ 
to  denote  that  there  are  continuous  projections  U,E  such  that  these  equalities  hold. 
The  assertion  that  &W  *  Z£  is  an  important  restriction  on  A  if  the  space  is 
infinite  dimensions. 

Lemma  1.1.  If  A  €  S/[X,Z) , >*tA) \  =  Xy,.5?(A)  =  ZE>  then  there  exists  a  right  inverse 
K  CMZg.Xj.u)  £f  A,  AK  =  I  on  Z£ , KA  *  I  -  U  on  X.  (see  Fig.  1.1) 


A  i?(A)«ZE 


^ ^ 

— i< 

>^(A)»XU 

Z.-E 

Figure  1.1 


Proof :  Since  A  is  one-to-one  from  Xjy  onto  Z£,  the  existence  of  K  is  clear.  The 
fact  that  K  is  bounded  follows  from  the  open  mapping  theorem. 

Suppose  A  is  a  linear  operator  as  in  Lemma  1.1,  A  is  a  Banach  space  denoting 
the  parameter  space  and  N  :  X  *  A  -*•  Z  is  a  C1  function  satisfying 

N(0,0)  ■  0,  DXN(0,0)  *  0  (1.1) 

We  want  to  discuss  the  solutions  of  the  equation 

Ax  =  N(x,X)  (1.2) 

for  (x,X)  in  a  neighborhood  of  (0.0). 

Using  the  projection  operator  E  in  Lemma  1.1,  we  can  rewrite  (1.2)  in  the 
equivalent  form 

EAx  «  EN(x,A),  (I-E)Ax  =  (I-E)N(x.X). 

If  we  let  x  »  y  ♦  z,  y  €  Xy,  z  €  Xjy,  and  use  the  fact  that  EA  *  A,(I-E)A=0,  Ax  =  Az 

and  K  is  a  right  inverse  of  A  on  Z  ,  we  obtain  the  equivalent  equations, 

t 

z  «  KEN(y+z,X)  (1.3a) 

0  ■  (I-E)N(y+z,X) .  (1.3b) 

One  can  use  the  Implicit  Function  Theorem  to  obtain  a  unique  solution  z*(y,x)  of 

(1.3a)  in  a  neighborhood  of  zero,  z*(0,0)  =  0,  Dyz*(0,0)  =  0.  For  x  =  y  ♦  z*(y,X) 

to  be  a  solution  of  the  equation  (1.2),  the  pair  (y,X)  must  satisfy 


G(y,X)  *  0 

G(y,X)  =  (I-E)N(y*z*(y,X),X). 


(1.4) 


The  function  G(y,X)  is  known  as  the  bifurcation  function.  The  above  procedure  for 
obtaining  solutions  of  (1.2)  is  an  application  of  the  alternative  method  and  is 
known  as  the  method  of  Liapunov-Schmidt  (LS  method) ■  It  is  summarized  in 


(1.2)  has  the  form  x  =  y  ♦  z*(y,X)  where  z*(y,X)  is  the  solution  of  (1.3a)  and  (v.X) 
satisfy  (1.4) . 

Several  specific  illustrations  of  Lemma  1.2  will  be  given  in  these  notes.  At 
the  same  time,  we  will  discuss  the  solutions  of  the  bifurcation  equation  and  relate 
the  analysis  to  dynamical  behavior. 

An  operator  A  :  X  ■*  1  with  closed  range  and  having  dim>f(A)  <  codim  4P(A)  <  •», 
is  called  a  Fredholm  operator  of  index  dim^A)  -  codim^P(A). 

In  the  applications.  Equation  (1.2)  often  arises  in  the  following  manner.  Sup¬ 
pose  M:  X  11  A  -*  Z  is  a  given  smooth  function  and  suppose  it  is  known  that  the  equa¬ 
tion  M(x, X)  *  0  has  a  solution  x  *  <p(X)  for  X  in  some  open  set.  One  can  study  the 
solutions  of  M(x,X)  *  0  near  4>(X)  by  letting  x  **•  ip(X)  +  x  to  obtain  a  new  function 

which  we  again  call  M  such  that  M(0,X)  =  0  for  X  in  an  open  set.  The  Taylor  series 

for  M  is  then 

M(x,X)  =  D(X)x  ♦  M(x, X) ,  M(x,X)  =  o(|x|)  as  |x|  -  0. 

If  the  operator  D(Xq)  has  a  bounded  inverse,  then  the  Implicit  Function  Theorem  im¬ 
plies  that  M(x,X)  =  0  has  a  unique  solution  x*(X)  in  a  neighborhood  of  (0,Xg), 
x*(XQ)  =  0.  Thus,  no  bifurcation  can  occur.  If  D(XQ)  does  not  have  an  inverse, 
then  there  is  the  possibility  of  bifurcation  near  (0,XQ).  In  this  case,  A  *  D(X0), 
N(x,X)  *  [D(X)-D(X0)]x  +  M(x,X). 

An  important  special  case  arises  when  X  is  a  scalar  parameter  and  D(X)  =  B  -  XC, 
where  B,C  are  bounded  linear  operators.  The  values  XQ  where  D(Xq)  is  singular  are 
then  eigenvalues  of  the  pair  of  operators  (B,C).  For  later  reference,  we  say  X^  is 

a  simple  eigenvalue  of  (B.C)  if  B  -  X^C  is  Fredholm  of  index  0  with  dimvfTB-  xoC)  *1  = 

codim i?tB-X0C)  and  OT(B-X0C)  «.S?tB-X0C)  =  Z. 

2.  Stable  and  unstable  manifolds.  In  this  section,  we  show  how  the  classical 
method  of  obtaining  stable  and  unstable  manifolds  for  an  hyperbolic  equilibrium 
point  is  a  special  case  of  the  LS  method  for  a  Fredholm  operator  with  A  :  X  -►  Z  with 
dim^tA)  <  <»  and  4?(A)  *  Z. 

Suppose  Aq  is  an  n  »  n  constant  matrix  with  Re  o(AQ)  i  0  where  o(A0)  is  the  spec¬ 
trum  of  Aq,  f  :  Jn  +  Rn  is  a  function  with  f(0)  *  0,  fx (0)  =  0.  In  order  to  con¬ 
struct  the  local  unstable  manifold  of 

(Ax)  (t)  «  f(x(t)),  (Ax)  (t)  =  i(t)  -  AQx(t),  (2.1)f 

we  consider  the -set  »  (initial  values  of  solutions  of  (2.1)  which  are  defined  and 
remain  in  a  "small"  neighborhood  of  zero  for  t  €  (-“,0]}.  The  local  stable  manifold 
Sj.  is  defined  similarly  on  (0,«).  Let  X°  *  (bounded,  uniformly  continuous  functions 
on  (-»,0|  to  Rn)  with  the  sup  topology  and  let  X*  *  (g  €  x"  :  g  €  X0}  with  the  C1 
topology.  Then  A  in  (3.1)f.  takes  x'  to  X(l  and  is  continuous  and  linear. 
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If  Un,S.  are  the  stable  and  unstable  manifolds  of  (2.1)0,  Rn  =  U.  ®  S  , 
u  u  10^^ 

7T  :  Fn  IL,  I  -  it  :  Rn  S  are  projections,  then  A  :  X  — ♦  X 
u  0  At 

-/f(A)  =  (g  €  X1  :  g(t)  =  e  0  *8(0),  t  1  0) . 

Thus,  dim/f^A)  <  One  can  define  a  projection  *  of  X1  onto  yf"(A)  by  the  relation 
(xrg)  (t)  =  e  °iTg(0),  t  ^  0.  It  is  a  classical  result  in  differential  equations  and 
easy  to  prove  that  ^A)  =  X°.  If  x  =  y  ♦  z,  y  *  e  €  xj  ,  2  6  xj  ~,  y(0)  =  £, 

then  the  LS  method  implies  there  is  a  unique  solution  xt(£,-)  =  y  +  zt(£,-)  €  X*  of 

1  *  1 

(2.1)^  in  a  sufficiently  small  neighborhood  of  zero  and  this  function  is  C  in  £, 
z*(0,-)  s  0,  D^z*(0,-)  =  0.  The  manifold  is  defined  by  (x  :  x  *  5  ♦  z*(£,0), 

£  small).  It  can  be  shown  that  any  solution  with  data  in  tl^  approaches  zero  as 
t  -*■  -■»  and  thus,  represents  the  local  unstable  manifold.  A  similar  argument  gives 
the  stable  manifold  S^.. 

Note  that  this  is  a  good  example  of  the  application  of  the  LS  method,  but 
there  is  no  bifurcation  equation  since  codim  J&(A)  =  0,  index  A  -  dim^ffA) . 

3.  Equilibrium  bifurcation.  In  this  section,  we  consider  the  (n+1 ) -dimensional 
system 

1  «  Cz  +  f (z , X )  (3.1) 

where  X  is  a  parameter  in  a  Banach  space  A,  f  :  Rn+1  x  A  -*Rn+1  is  continuous 
together  with  derivatives  up  through  order  k  >  1, 

f(0,0)  =  0,  3f(0,0)/3z  =  0  (3.2) 


(3.3) 


and  B  is  an  n  x  n  matrix  with  Re  ct(B)  t  0,  where  o(B)  is  the  spectrum  of  B.  If 
z  *  (x,y),  f  -  (g,h),  with  x  EH,  g  E  R,  Eq.  (1.1)  can  be  written  as 


x  -  g(x,y,X) 
y  »  By  ♦  h(x,y,X) . 


(3.4) 


Our  objective  is  to  discuss  the  bifurcation  and  stability  of  equilibrium 
solutions  of  Eq.  (3.1)  in  a  neighborhood  of  (z,X)  *  (0,0).  The  equilibrium  solutions 
are  easily  obtained  by  the  LS  method  applied  to  the  equation  Cz  +  f (z,X)  »  0.  This 
is  equivalent  to  applying  the  Implicit  Function  Theorem  to  obtain  a  solution  <p(x,X) 
of  the  equation 


Bxp  +  h(x,ip,X)  =  0 


(3.51 


in  a  neighborhood  of  (x,X)  =  (0,0)  which  satisfies  <p(0,0)  *  0,  S«i{0,0)/3x  =  0.  The 
equilibrium  points  of  Eq.  (3.1)  in  a  neighborhood  of  (z,X)  =  (0,0)  are  then  given  by 
(x,<p(x,X))  whdre  (x,X)  satisfy  the  bifurcation  equation 
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G(x,X)  «  0 

G(x,X)  «  g(x,q>(x,X),X). 


(3.6) 


We  remark  that  the  function  <p(x,X)  depends  only  upon  the  vector  field  h(x,y,\)  [not 
on  g(x,y,X)]  and  has  the  same  smoothness  properties  as  h(x,y,X). 

The  equilibrium  points  of  Eq.  (3.1)  can  also  be  expressed  in  terms  of  the 
center  manifold.  In  fact,  there  is  a  center  manifold  =  {(x,y)  :  y  =  ii<(x,X),  x 
in  a  neighborhood  of  zero)  for  X  in  a  neighborhood  of  zero.  The  flow  on  is 
given  by  (x(t),y(t))  =  (x(t)  ,i|f(x(t)  ,X))  where  x(t)  satisfies  the  equation 


x  =  G(x,X) 

C(x.X)  =  g(x,<Kx.X).X). 


(3.7) 


The  equilibrium  points  of  Eq.  (2.1)  in  a  neighborhood  of  (z,X)  =  (0,0)  are  given  by 
(x.iKx.X))  where  G(x,X)  «  0. 

We  remark  that  the  function  <i(x,X)  depends  on  both  of  the  vector  fields 

g(x,y,X),  h(x,y,x)  and,  generally,  is  not  as  smooth  as  g  and  h.  In  particular, 

00  00 

i|/  will  generally  not  be  C  or  analytic  even  when  g,h  are  C  or  analytic. 

Since  a  center  manifold  has  a  hyperbolic  structure  (each  point  on  looks 
like  a  saddle  point),  the  flow  on  the  center  manifold  gives  a  complete  description 
of  the  flow  defined  by  Eq.  (3.1)  in  a  neighborhood  of  (z,X)  =  (0,0).  Even  though 
the  bifurcation  function  G(x,X)  in  (3.6)  was  constructed  without  mentioning  dynamical 
behavior,  it  is  an  interesting  fact  (stated  precisely  in  Theorem  3.1  below)  that  the 
flow  defined  by  the  equation 

x  =  G(x,X)  (3.8) 

is  equivalent  to  the  flow  defined  by  Eq.  (3.7).  Thus,  the  qualitative  properties  of 
the  flow  can  be  determined  without  knowing  the  center  manifold.  The  advantage  of 
this  observation  lies  in  the  fact  that  G(x,X)  is  as  smooth  as  the  original  vector 
field  and  is  easier  to  calculate  approximately.  Up  to  this  point,  the  parameter  X 
has  played  no  role  and  it  plays  no  role  in  the  following  theorem.  However,  it  will 
be  used  in  a  significant  way  in  the  applications  of  Theorem  3.1. 

Theorem  3.1.  The  vector  fields  G(-,X),  G(-,X)  are  equivalent  in  a  neighborhood  of 
zero  for  X  in  a  neighborhood  of  zero;  that  is,  there  is  a  homeomorphism  of  a  neigh¬ 
borhood  of  x  »  0  mapping  the  orbits  of  (3.7)  onto  the  orbits  of  (3.8)  preserving  the 
s ense  of  di rection  in  time. 


Proof:  The  functions  G(x,X),  G(x,X)  must  have  the  same  set  of  zeros  in  a  neighbor¬ 
hood  of  (x,X)  =  (0,0).  The  essential  element  of  the  proof  of  the  theorem  and  the 
only  part  that  will  be  given  is  to  show  hat  G(x,>),  G(x,X)  have  the  same  sign 


between  zeros. 


Suppose  this  is  not  the  case;  that  is,  there  is  an  xQ  such  that 


C(x0,A)  C(x0,A)  <  0. 


By  making  a  small  perturbation  of  the  vector  f(-,A)  in  Eq. 


(3.1),  one  obtains  new  functions  G^(x,A),  Gj(x,A)  such  that  Gj(x0>A)  Gj(xq,A)  <  0 


and  the  first  zero  Xj  before  x^  is  simple.  Now  make  another  perturbation  of  f (•,'.) 
in  Eq.  (3.1)  by  replacing  g  by  g  +  c  with  e  small  >  0.  The  new  bifurcation 
function  G,(x,A,e)  =  Gj(x,A)  +  e.  Also,  using  the  theory  of  center  manifolds,  one 
can  show  that  the  vector  field  on  the  center  manifold  G, (x,A,c)  -  Gj(x,A),  + 
o(x,A,e)e  where  o(x,A,e)  >  0.  Since  Xj  is  a  simple  zero  of  Gj(x,A),  Gj(x,A),  this 
implies  there  are  functions  x2(X,e),  x2(A,e),  x2(A,0)  =  ^  »  x2(X,0),  such  that 
G2(x2(X,e)  ,A,e)  =  0,  G2(i2(A,e),A,e)  =  0  and  such  that  [dx2(A,0)/de]  •  (dx^  (>  ,0)./dE] 

<  0,  which  implies  x2(A,e)  jt  x2  (A ,  c),  e  >  0.  But  this  is  a  contradiction  since  the 
bifurcation  function  and  the  vector  field  on  the  center  manifold  must  always  have 
the  same  zeros.  This  implies  G(x,A)  and  G(x,A)  have  the  same  sign  between  zeros. 


As  a  first  application,  we  give  a  result  on  bifurcation  from  a  simple  eigen¬ 


value  for  one  parameter  families  of  vector  fields, 
write  f (z , A)  as 


Suppose  A  £k,  f(z,A)  is  C  and 


f (z,A)  =  C^AJz  ♦  F(z,A) 


(3.9) 


F (0, A)  =  0,  3F(0,A)/3z  =  0  for  all  A 
where  Cj(A)  is  an  (n+1)  x  (n+1)  matrix,  Cj(0)  =  0.  Suppose  that 


A  0 

:.(*)  = 

1  Lo  B.(A)J 


(3.10) 


and  Bj(A)  is  an  n  *  n  matrix.  The  bifurcation  function  G(x,A)  for  this  special 
case  has  the  form 

C(x,A)  =  a(A)x  +  G(x,  A)  ,  (3.11) 

where  G(x,A)  *  0(|x|2)  as  x  +  0,  a(0)  *  0,  a  (0)  =  1  from  (3.9),  (3.10).  The 
Implicit  Function  Theorem  and  Theorem  3.1  imply 

2 

Theorem  3.2.  If  A  €  R,  f  is  C  and  satisfies  (3.9),  (3.10),  then  there  is  a  neigh¬ 
borhood  V  of  z  =  0,  W  of  A  =  0,  U  of  x  =  0  and  a  C1 -function  A*  ;  U  ■+•  R  such  that 
A*(0)  =  0  and,  for  A  €  W,  the  equilibrium  solutions  of  (3.1)  in  V  are  given  by 
{0}  U  {(x,q>(x,A))  :  A*(x)  =  A).  A  solution  (x,ip(x,A))  is  asymptotically  stable 
(unstable)  on  the  center  manifold  if  (x-x)G(x,A)  <  0  (>0)  near  x. 

The  stability  properties  are  easy  to  determine  from  the  curve  r  :  A  =  A*(x), 
x  €  U.  In  fact,  since  3G(0,A)/3x  =  A  +  o(|x|)  as  A  ■*  0,  it  follows  that,  for  a 
fixed  A  i  0,  sign  xG(x,A)  =  sign  A  for  x  near  zero.  This  implies  that  the  sign  of 
G  is  as  indicated  in  Figure  3.1.  This  yields  stability  of  all  equilibrium  points 
immediately  with  the  situation  as  shown  in  Figure  (3.1).  At  the  bifurcation  points 
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there  is  always  a  transfer  of  stability  on  the  center  manifold. 

i 


Figure  3.1 

2 

As  another  illustration  suppose  f  in  (3.1)  is  C  and  X  is  a  parameter  in  A 
and  there  is  a  6  t  0  such  that 

G(x,0)  *  Bx2  ♦  o ( | x | 2)  as  x  ■*  0.  (3.12) 

Since  (p(0,0)  *  0,  3<p(0,0)/3x  =  0,  it  follows  tha  :  6  is  given  explicitly  by 

28  *  32g(0,0,0)/3x2.  The  Implicit  Function  Ther rem  implies  there  is  a  unique  C1- 

function  a(X),  a(0)  =  0,  such  that  3G(a  (X) ,  X)/3x  =  0  for  X  small.  Then 

G(a(X)  +  x,X)  =  yQ(X)  ♦  y2(X)x2  ♦  °(|x|2)  as  |x|  ♦  0 

where  Yq(0)  =  0  and  y^(0)  =  8  j*  0.  Theorem  3.1  therefore  implies,  near  (z,X)  ■ 
(0,0),  that 

(i)  yq(A)  >  0  =*>  no  equilibrium  solution  of  (3.1). 

(ii)  Yfl(X)  =  0  one  equilibrium  solution  of  (3.1)  stable  from  one  side 

and  unstable  from  the  other  side  on  the  center  manifold. 

(iii)  Y0(i)  <  0  implies  two  equilibrium  solutions  of  (3.1),  hyperbolic 
with  one  stable  and  the  other  unstable  on  the  center  manifold. 

The  surface  in  A  defined  by  T  *  {'*  :  yQ(X)  =  0}  is  called  the  bifurcation  surface. 

This  bifurcation  is  referred  to  as  the  saddle-node  bifurcation.  The  name  comes  from 

the  fact  that,  for  z  in  R  ,  a  saddle  and  node  coalesce  and  disappear  as  r  is 
crossed. 

As  another  illustration,  suppose  the  vector  field  in  (3.1)  is  C3  and  there  is 
a  8  i  0  such  that 

G(x,C)  =  Bx3  ♦  o(|x|3)  as  |x|  ■*  0.  (3.13) 

2  2 

Then  there  is  a  unique  point  a(X),  a(0)  =  0,  where  3  G(x,X)/3x  =0  in  a  neighbor¬ 

hood  of  zero.  x  •  a(X)  ♦  x,  then 

G(a(X)  ♦  I.X)  =  Y0(X)  ♦  Yj (X)X  +  Yj(X)I3  ♦  G(x,X)  (3.14) 

where  G(x,X)  =  o(|x|3)  asx-*0,  Yq(0)  =  Yj  (0)  *  0,  Yj(0)  =  8  i  0. 

It  is  no  loss  in  generality  to  assume  Yj(X)  =  1  in  discussing  the  solutions  of 

(3.14)  in  a  neighborhood  of  (x,X)  =  (0,0).  F.ven  in  the  case  when  C.(x,X)  -  0  in 

(3.14) ,  a  complete  discussion  of  the  zeros  of  the  cubic  requires  two  parameters. 
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This  corresponds  to  a  codimension  two  bifurcation.  A  generic  two  parameter  family 
of  parameters  A  *  (A^.A^)  should  satisfy  det  [3 (YQ,y j)/3  Pq.A  j )  ]  t  0  at  X  =  0.  If 
we  assume  this  is  the  case,  then  we  can  introduce  new  parameters  u  =  (u0>.jj)  instead 
of  A  to  obtain  the  equation 


H(X,ll)  *  ♦  DjX  +  x3  +  H(x,u)  =  0 


(3. IS) 


where  H(x,u)  =  o(|x|3)  as  x  -*  0.  The  bifurcation  curves  in  u-space  correspond  to 

multiple  solutions  of  (3. IS);  that  is,  H(x,w)  =  0,  3H(x,u)/3x  =0.  In  a  neighbor- 

—  —2  |  — ,  2 

hood  of  zero,  this  determines  u.,u  as  functions  of  x,  u  =  -3x  ♦  o(|x|  ), 

“3  _  __  u  l  _  * 

u  =  2x  +  o ( | x |  )  as  x  •*  0.  Eliminating  x  from  these  equations,  one  obtains  a 

0  3  2 

cusp  in  u-space  which  is  approximately  given  by  u^  *  -(27/4)u0-  The  bifurcation 

diagrams  for  Equation  (3.15)  are  shown  in  Fig.  3.2  with  the  number  of  solutions  as 

indicated.  ....  i  ,  u~>0 


Figure  3.2 

For  some  problems,  one  may  have  a  parameter  A  €  A  for  which  rank  l3 (Yq> Y j ) /3A] 
<  2  at  A  *  0.  In  this  case,  the  bifurcation  diagram  will  be  different  since  the  set 
{ (Yq (A ) , Y j (A) ) , A  6  A)  will  be  a  curve  in  (u0>Uj)-space.  The  number  of  solutions  and 
the  bifurcations  for  Equation  (3.14)  will  be  determined  by  how  this  curve  intersects 
the  bifurcation  curve  in  Fig.  3.2  as  A  varies. 

If  G(x,0)  =  Bxq  *  o(|x|q)  as  |x|  -►  0,  where  8^0,  then  the  two  previous  ex¬ 
amples  show  that  only  one  parameter  is  needed  to  describe  all  bifurcations  when 
q  «  2  and  two  are  needed  to  describe  the  bifurcations  when  q  =  3.  For  q  >  3,  one 
can  show  that  q  parameters  are  needed.  The  bifurcation  surfaces  are  ex.remely 
complicated  and  will  not  be  discussed  here. 

Theorem  3.1  is  valid  for  several  types  of  partial  differential  equations  and 
functional  differential  equations.  As  is  often  the  case,  the  theoretical  results 
are  very  simple  to  state,  but  there  are  many  technical  problems  in  making  the 
applications.  We  briefly  indicate  how  this  can  be  done  for  retarded  functional 
differential  equations  and  for  abstract  evolutionary  equations  which  include 
parabolic  equations. 

Let  A  be  the  infinitesimal  generator  of  an  analytic  semigroup  T(t),  t  >  0, 
on  a  Banach  space  X.  Choose  b  real  so  that  Re  a  (A)  ■>  0,  A  »  A  ♦  bl.  One  can 
define  fractional  powers  A°  of  A  for  any  a  £  R.  For  a  <_  0,  a  is  one-to-one  and 
a  bounded  linear  operator  on  X.  For  a  ^  0,  consider  the  set  X  c  X  consisting 


of  all  x  €  X  for  which  lA  xj  < 


For  x  £  X ,  define  |  x  [ 


1 1  *  | A  ’ x  j .  Then 


10 


|*|  is  a  norm  on  Xa  and  Xa  is  a  Banach  space  with  this  norm.  Also,  the  inclusion 
map  taking  Xa  into  X  is  continuous.  Now  suppose  that  f  :  Xa  x  X  is  a  smooth 
function  for  some  0  <_  a  <  1  and  consider  the  evolutionary  equation 

u  +  Au  =  f(u,X) .  (3. 16) 

Suppose  that  Eq.  (3.16)  generates  a  strongly  continuous  semigroup  T^(t),  t  >_  0,  on 
X  .  The  choice  of  a  depends  upon  the  specific  function  f. 

Suppose  0  €  o(A)  is  a  simple  eigenvalue  of  A  and  there  is  a  6  >  0  such  that 
|Re[o(A)  n  (0}]|  >_  5.  Then  the  space  X°  can  be  decomposed  as  a  direct  sum  Xa  =  V  •  h' 
with  V,N  invariant  under  T\(t)  where  V  is  one  dimensional  and  spanned  by  a  unit 
vector  Up  satisfying  AuQ  *  0.  If  u  ■  v  *  w,  v  £  V,  w  £  K,  f  *  g  ♦  h,  g  £  V,  h  €  W, 
then  (3.16)  is  equivalent  to 

v  *  g(v*w,X) 

(3.17) 

w  »  -Aw  ♦  h(v»w,X) . 

The  operator  A  restricted  to  N  has  a  bounded  inverse  which  we  write  as  A-1 
taking  X  into  X.  Since  D(A)  c  D(A°)  fur  0  ^  a  <  1,  it  follows  that  the  equation 

-A*P  ♦  h(v*ip,  X)  «  0,  ip  £  K,  (3.18) 

is  equivalent  to  the  equation 

<P  ♦  A*lh(v««P,k)  -  0,  < PEW.  (3.19) 

Consider  the  operator  A  1  as  mapping  X  into  Xa,  0  <_  a  <  1.  Then  A1  is  a  bounded 
linear  operator  from  W  to  X°  since  A"1+a  is  a  bounded  linear  operator  on  X  (-l+a<0) 
and  |A_1u|  *  j  A- 1  u |  ♦  |  Ac*A_1u|  ♦  |  A- 1  u |  ♦  | A- 1  ( I *bA_  1 ) u  |  for  any  u  €  D(Aa)  fl  W. 

The  Implicit  Function  Theorem  implies  there  is  a  unique  function  ip(v,X), 
v(0,0)  *  0,  satisfying  (3.19).  The  equilibrium  points  of  (3.16)  are  therefore  given 
by  u  «  v  ♦  ip(v,X) ,  where  u,X  must  satisfy  the  bifurcation  equation 

G(v,X)  *  0,  G(v,X)  ■  g  (v+ip(v , X)  , X)  . 


The  flow  on  the  center  manifold  of  (3.16)  is  equivalent  to  the  flow  defined  by  the 
scalar  ordinary  differential  equation  v  =  G(v,X). 

As  a  specific  example,  consider  the  equation 


u=u  ♦  jiu  -  u  ,  0  <  x  <  n 

t  XX  ’ 

u  =  0  at  x  =  0,*. 


(5.20) 


Let  X  *  L2(0,*),  *l’2  *  (<P  €  X  :  <P  €  X,  <P(0)  »  <P(*)  *  0)  and  define  Aip  *  -<P  with 

1  2  2  2°  x  xx 
D(A)  *  N  ’  fl  «  ’  .  One  can  show  that  the  operator  A  generates  an  analytic  semi- 

^  l  /2  1  1  1 / 1 

group  on  X.  It  is  also  possible  to  show  that  D(A  )  *  N  '*  *  X  ,  the  operator  A 

I)  22 

has  compact  resolvent  and  o(A)  consists  of  simple  eigenvalues  (1,2  . n  ,...)  with 


1 1 


the  eigenfunction  (2/i’)1'/2sin  nx  corresponding  to  n2.  If  F(ip,p)(x)  =  uip(x)  -  ip3(x), 

0  <  x  <  n,  then  one  can  show  that  F  :  X1^2  *  F  X  is  an  analytic  function  of  <p,p. 

—  —  i/2 

This  is  enough  to  conclude  that,  for  any  tp  €  X  ,  there  is  a  function  u(t,x,ip), 

1/2 

defined  on  some  interval  0  ^  t  <  t  ,  u(t,-,<p)  €  X  ,  0  <_  t  <  t^,  u(0,-,<P)  *  <P, 
u(t,x,<p)  is  continuous  in  (t,x,<p)  and  satisfies  (3.20)  for  0  <_  t  <  t  . 


For  Equation  (3.20),  there  is  a  Liapunov  function 


V(cp) 


■fij 

J0 


l^(x)  -  XH(cp(x) ) ]dx 


for  <p  € 


X1/2  and  H(u)  =  pu2/2  -  u4/4. 


In  fact, 


V (u(t ,  • ))  =  £  V (u(t ,  • ) ) 


u2(t,x)dx  <_  0. 
0  t 


One  can  use  this  inequality  to  show  that  u(t,x,<p)  is  defined  for  t  ^  0;  that  is, 

t  =  ».  If  one  now  defines  T  (t)ip  =  u(t,-,ip),  then  T  (t),  t  >  0,  is  a  strongly 
<p  l/2  v  v 

continuous  semigroup  on  X 

Let  us  now  determine  the  nature  of  the  bifurcation  near  the  simple  eigenvalue 
«  n2.  Let  p  =  n2  ♦  X  and,  for  u  €  X^2, 

u(x)  «  (2)1/2v  sin  nx  +  w(x),  vdif(-)1/2  f  u(x)sin  nx  dx 

n  "  J0 

f  w(x)sin  nx  dx  -  0. 

JO 

If  Uu  =  v  sin  nx,  then  U  is  a  projection  onto  ^"(A-p  I)  and  (I-U)is  a  projection 
onto  .SffA-p^I).  With  this  decomposition,  the  equilibrium  solutions  of  (3.20),  that 
is  ut  «  0,  must  satisfy  the  equations 

1/2  3 

Xv  sin  nx-  U[(2/n)  '  v  sin  n  •  +  w(-)]  -  0 

w  +  *  Jw  "  (I-^H(2/r)^2  sin  n  •  *  w(-)]3  «  0, 

with  w  »  0  at  x  -  0,ir.  This  last  equation  has  the  solution  w(v,X)  £  X  satisfying 
w(0,0)  =  0,  3w(0,0)/3v  »  0,  I  w(v,X) (x)sin  nx  dx  =  0.  Therefore,  the  equilibrium 
solutions  of  (3.20)  near  (u.p?  *  0,0)  are  given  by 

u(x)  -  (2)1/2  v  sin  nx  ♦  w(v,X)(x) 

TT 

with  (v,X)  satisfying  the  bifurcation  equation 
G(v,X)  =  0 

G(v,X)  *  Xv  -  (-)1/,2f  (sin  nx)[(-)I//2  v  sin  nx  +  w (v ,X)  (x) ] 3dx. 

"JO  ” 

1  t  I 

The  function  G(v,X)  *  y0(X)v  -  Yj(X)v  ♦  o(|v|  )  as  v  -»  0  where  -y0(0)  s  1, 

Vj(0)  *  3/tt2.  Thus,  the  bifurcation  curve  x  *  X*(v)  satisfying  G(v,X*(v))  !  0  is 
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2  2  2 

given  by  X  =  (3/r  )v  +  o(v  )  as  v  -»  0  (see  Fig.  3.3).  At  each  bifurcation,  there 


vl 

1 


Figure  3.3 


is  a  transfer  of  stability  on  the  center  manifold.  However,  for  the  complete  flow, 
the  only  stable  equilibrium  points  are  the  ones  that  occur  at  the  first  bifurcation. 
This  is  because  the  zero  solution  has  an  unstable  manifold  of  dimension  >  1  at  the 
other  points. 

As  another  illustration,  we  consider  a  retarded  functional  differential  equa¬ 
tion.  Suppose  r  >_  0  is  a  given  constant,  C  =  C([-r,0],  Rn) ,  L  :  C  +  Rn  is  a  con¬ 
tinuous  linear  operator,  f  :  C  x  A  •*  Rn  is  a  given  C^-function,  f(0,0)  =  0, 

D^f(<p,X)  =  0  when  (<p,X)  =  (0,0).  For  a  given  function  x  :  f-r,A)  -►  Rn  and  a  fixed 
t  €  [0,A),  we  let  x^  designate  the  function  from  [-r,0]  to  Rn  given  by  xt(e)  = 
x (t+6 ) ,  -r  e  £  0,  and  consider  the  retarded  functional  differential  equation 

i(t)  -  Lxt  ♦  f(xt,X).  (3.21) 

If  v>  €  C,  let  x((p)(t)  be  the  solution  of  (3.21)  with  xQ(ip)  *  ip.  We  suppose  all 
solutions  are  defined  on  [-r,«)  and  define  T^(t)<p  =  xt(ip).  Then  T^(t),  t  ^  0,  is  a 
strongly  continuous  semigroup. 

An  equilibrium  point  for  (3.21)  is  a  solution  x  defined  for  all  t  ED,  xQ  - 
and  T^(t)i^  *  i 1)  for  all  t  €  R.  This  is  equivalent  to  x(t+9,i(i)  =  il(0)  for  all  t  6  R. 
Since  x (t , ip )  is  continuously  differentiable  for  all  t  ER  (this  is  a  consequence  of 
the  fact  that  it  exists  for  t  €  R) ,  this  implies  x(t)  =  0  for  all  t  £R;  that  is, 
x(t)  =  b,  a  constant,  for  all  t  6  R  and  iji(e)  =  b,  -r  <_  e  <_  0.  Therefore,  the 
equilibrium  points  of  (3.21)  are  constant  functions  \ji  such  that 


Lip  ♦  f  (ip ,  X )  =  0 

(3.22) 

i(e)  =  o,  -r  <_  e  <_  o. 

The  linear  equation 

x(t)  *  L(xt)  (3.23) 

also  generates  a  strongly  continuous  linear  semigroup  S(t),  t  >_  0,  on  C.  It  can  be 
shown  that  the  infinitesimal  generator  A  of  S(t)  is  given  by  (Atp)  (e )  =  ip(e), 

"r  1  9  L  0,  with  °(A)  =  (»  £  C  :<p0  €  C,  ip(0)  =  L(<p)}.  Also,  o(A)  consists  only  of 
point  spectrum,  and,  if  L(<p)  «  (dn (0)]ip(e) ,  then  u  €  o(A)  if  and  only  if 

*r  f° 

det  d(p)  •  0,  fi(y)  -  ul  -  eu0dn(6).  (3.24) 


i 
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Also,  if  m  €  o(A),  then  ./f'(A-uI)  is  the  span  of  the  functions  e  b,  -r  £  6  <_  0, 
where  A(u)b  *  0. 

Suppose  0  £  o(A)  is  simple  and  Re  u  t  0  for  v  €  o(A)  v  (0).  Then  it  follows 
immediately  from  (3.23)  that  there  is  a  6  >  0  such  that  |Re  w|  >_  6  for  v  e  ofA)"v(o> 
Our  first  objective  is  to  determine  the  equilibrium  points  of  (3.21)  near  x  =  0  for 
\  near  zero. 

To  characterize  the  equilibrium  solutions  of  (3.21)  as  the  solutions  of  an 
operator  equation,  we  need  to  enlarge  the  space  of  functions  that  are  being  con¬ 
sidered.  let  XQ  be  the  n  x  n  matrix  function  on  [-r,0]  defined  by  XQ(e)  =  0, 

-r  <_  e  <  0,  X0(0)  *  I,  the  identity,  and  define  PC  =  C  ®  (Xp)  where  (XQ)  denotes 
the  space  of  XQ.  A  function  €  PC  is  represented  by  a  pair  (<p,b)  where  <p  €  C  and 
b  €  R  ;  the  function  V  is  uniformly  continuous  on  [-r,0),  i|«(e)  =  ip(6),  6  6  [-r,0), 
and  >(((0)  *  ii>(0)  ♦  b;  that  is,  <(i  has  a  jump  of  b  at  0.  Consider  A  as  a  map  from 
C  to  PC  with  D(A)  *  (<p  €  C  :  <9  €  C}  and 

(Alp)  (6)  =  <p(9),  -r  i  6  <  0 

(Alp)  (0)  *  UP. 

It  is  now  easy  to  see  from  (3.22)  that  ♦  is  an  equilibrium  point  of  (3.21)  if  and 
only  if 

A*  ♦  X0fOM)  *  o  (3.25) 

We  now  apply  the  LS  method  to  Eq.  (3.25)  under  the  assumption  that  0  is  a 
simple  eigenvalue  of  A.  We  know  that  >^(A)  =  (q>  €  C  :  >P  =  fkPg,  VPQ (8 j  =  a,  -r  <_  6  50, 
where  a  is  a  constant  n-vector,  A (Oja  =0).  To  characterize  the  range  of  A,  we 
solve  the  equation  Atp  *  1)/  for  €  PC.  This  relation  is  equivalent  to  tp(0)  =  <1(6), 

-r  <  0,  and  UP  *  *(0);  that  is, 

f® 

ip(6)  •  ip(0)  ♦  *(s)ds 

Jo  0  e 

M0)ip(0)  *  *(0)  -  f  [dn(9)]|  i|/(s)ds. 

J-r  ;0 


Let  b 


be  any  nonzero  row  vector  such  that  bi(0)  * 

1  0,  o.0(6)  •  b, 

(»„.*)  ■  bK>(0)  -  f  b[dn(6)][ 
u  J-r  J0 


0.  Then  < l>  C^KA)  if  and  only  if 
-r  1  0  1  0 

(3.26) 


iKsjds. 


Choose  b  so  that  (ao>%)  *  1  with  in  (3.26)  and  iPQ  a  basis  for  ^(A) .  The 
space  PC  can  be  decomposed  as 


PC  *  ^A)  •  #(A) 

O'  €  pc,  *  •  ftP0  ♦  n,  8  »  (a0,*)  f  R.C^p.n)  «  0. 


(3.27) 
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These  same  formulas  also  give  a  decomposition  of  C  as  C  =  P  ®  Q,  P  =  (ip=8q>0,8€R}t 
Q  =  {«P  :  (a0'<P)  *  0). 

The  computations  above  show  that  A  *  exists  and  is  bounded  as  a  map  from  4P(A) 
into  Q.  We  can  now  apply  the  LS  method. 

If  Xq  «  iPpd,  <5  =  (a0»x0)>  Xq  '  xo  "  ^o5,  then  that  t)>  =  6q>0  +n, 

(aQ,n)  «  0,  is  an  equilibrium  point  of  (3.21)  if  and  only  if 

bf((PQ6+n,X)  =  0 

n  +  A^xjjfOPjjB+n.A)  =  o. 

The  last  equation  has  a  unique  solution  n(B,X)  in  a  neighborhood  of  (B,X)  =  (0,0), 
n(0,0)  «  0,  3n(0,0)/38  =  0.  The  equilibrium  solutions  of  (3.21)  are  given  by 
i)/  »  <pQ8  ♦  n(6,X)  where  (8,X)  satisfy  the  bifurcation  equation 

G(B.X)  *  bf(«P0B*n(8,X),X)  «  0.  (3.28) 

One  can  also  show  that  the  flow  on  the  center  manifold  of  (3.21)  is  equivalent  to 
the  flow  defined  by  6  =  G(6,X). 

As  an  example,  consider  the  linear  equation 

x(t)  «  -x(t)  ♦  x(t-l) 

which  has  g(X)  »  X-l-exp(-X).  2ero  is  a  simple  eigenvalue  and  (aQ,\(i)  * 

b(®(0)  ♦  |  ®(s)ds] .  If  we  choose  ^(9)  «  1,  -r  6  <.  0,  then  (ctg,<Pg)  =  1  if 

b  =»  1/2.  If  we  choose  f(ip,X)  =  X  +  <p2(0),  then  the  bifurcation  equation  (3.28) 

2 

is  equivalent  to  X  -  (6+n (8, X) (0) )  =  0.  Thus,  there  will  be  two  equilibrium  points 
for  X  >  0  and  no  equilibrium  points  for  X  <  o.  The  same  remark  is  true  if  f(<P,X)  * 

X  -  «P2(-1). 


Bifurcation  from  two  purely  imaginary  roots.  In  this  section,  we  consider  the 


equation 


z  *  Cz  ♦  Z(z,y) 


where  u  1  E  is  a  parameter  in  a  Banach  space  E,  z  £Rn+2,  Z  is  a  C^-function, 
k  ^  2,  C  is  an  (n*2)  *  (n*2)  matrix  satisfying 


Z(0,u)  =  0,  3Z(0,0)/3z  =  0,  for  all  y. 


[A 

°1. 

A  « 

r  °  n 

to 

Bj 

l-i  oj 

with  the  eigenvalues  of  B  having  nonzero  real  parts  (for  much  of  the  discussion 
2no 

(e  -I)  being  nonsingular  is  sufficient). 


Our  objective  is  to  discuss  the  periodic  orbits  of  (4.1)  in  a  neighborhood  of 
(z,u)  «  (0,0).  These  orbits  must  have  period  approximately  2"  and  be  close  to  the 
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2ir-periodic  solutions  of  the  linear  equation  z  =  Cz.  If  t  i — >  Bt,  where  2*/B  is  the 
period  of  the  solutions  to  be  determined,  then  S  should  be  close  to  1  and  the  new 
equation  for  z  is 

z  =  6Cz  +  BZ(z,y).  (4.3) 

We  will  apply  the  LS  method  for  2»-periodic  solutions  of  this  equation.  If 
z  *  (x,y),  Z*(HY)  ,  x,X  £R2,  then  (4.3)  is  equivalent  to 

x  -  Ax  =  (B-l)Ax  ♦  BX(x,y,w) 

(4.4) 

y  -  By  =  (S-l)By  +  BY(x,y,u). 

Let  pd  be  the  space  of  2n-periodic  (n+2-vector  functions  which  are  continuous  with 

1  0 
the  sup  topology  and  p7  be  those  2ir-periodic  (n+2) -vector  functions  in  p.  whose 

first  derivatives  are  also  in  p^with  the  C  topology.  If  Si/:  p7Tt  ■*  p  is  defined 

on  functions  f  =  (g,h)  by 

a/  f  =  (g-Ag,h-Bh) 

thenjj^is  a  continuous,  linear  operator  and  ^<e/)  *  {(eAtb,0),  0  €  Rn,  b  £  R2).  The 
classical  Fredholm  alternative  for  linear  periodic  systems^implies  &<s/)  -  (f=(g,h)  : 

[  e  Atg(t)dt  *  0).  Define  E  =  diag(E,I),  Eg  =  eAtD  *[  e  Atg(t)dt,  D  *  2*1. 

J0  0  0 
Then  gHjfi  =  (I-E)p"w. 

One  can  now  apply  the  LS  method.  Before  doing  this,  it  is  convenient  to  ob¬ 
serve  that  -*'s/j  is  determined  from  a  phase  shift  on  the  one  dimensional  subspace 
a(cos  t,-sin  t,0),  a  €  R.  Since  (4.4)  is  autonomous,  we  may  therefore  fix  the 
element  in  yf'TA)  as  a(cos  t,-sin  t,0).  An  application  of  the  LS  method  yields  a 
z*(8,a,y)  €  p2n  of  the  form 

z*  =  z*(B,a,w)  «  (a«Hx*(B,a,y),y»(B,a,y)) 

f  -At  (4-5) 

<p(t)  *  (cos  t,-sin  t),  e  x*(6,a,v)  (t)dt  -  0 

;0 

and  x*  «  x*(B,a,u),  y*  =  y*(8,a,u)  satisfying  x*(B,0,y)  =  0,  y*;8,0,u)  =  0, 

3x*(l ,0,0)/3a  =  0,  3y* (1 ,0,0)/3a  *  0  and  the  equations 

x  -  Ax  =  (I-E)  [  (8-l)A(aip*x)  +  BX(a4>»x,y,u) 


y  -  By  =  (8-1) By  ♦  6Y(a<JHx,y,u) . 

The  function  z*  will  satisfy  (4.1)  if  and  only  if  E[(6-l)A(3<P+x*)  ♦  6X(aip*x*,y*,u)]=0. 
This  latter  equation  is  equivalent  to 

G(8,a,u)dff[  e~At[  (8-l)A(a(p(t)+x*  (t )  )*BX  (a<b(t  )+x*  (t )  ,y*(t)  ,y))dt  =  0.  (4.7) 

'0 

This  represents  two  equations  in  the  parameters  (8,a,u).  However,  it  is  easy  to  see 


that  one  of  these  equations  can  always  be  solved  for  8  (this  determines  the  period) 
as  a  function  of  a,y.  In  fact  a  simple  calculation  shows  that  e'AtAq>(t)  =  (0,-1). 
Thus,  if  G  «  (gi*C2^  in  (4-7)»  then  G2(8.a,y)  satisfies  G2(8.0,u)  =  0,  3G2(1,0,0)/ 
3a  *  0  and  3^G2 (1 ,0,0)/3a38  »  -1.  This  implies  that  the  function  G2(8,a,y)/a  has 
a  unique  solution  8  *  8*(a,u)  in  a  neighborhood  of  (e,a,y)  =  (1,0,0),  8*(0,0)  =  1. 
Thus,  equations  (4.7)  are  equivalent  to  the  scalar  equation 


G(a,y)  =  0 


(4.8) 


G(a,y)  =  Gj  (8* (a,y) ,a,y) . 

We  refer  to  the  function  G(a,y)  as  the  bifurcation  function  for  periodic  orbits  of 
(4.1)  near  (i,y)  «  (0,0). 


Summarizing  the  above  remarks,  we  see  there  is  a  periodic  orbit  of  (4.1)  near 
(z,  y)  *  (0,0)  if  and  only  if 


z  »  (a  ♦x*(8*(a.y),a,w),y*(8*(a,y),a,u))  (4.9) 

where  x*(8,a,y),y*(8,a,u)  satisfy  (4.6),  and  (a,y)  satisfy  (4.8).  The  period  of 
z  is  2ir/8*(a,u).  All  functions  can  be  approximated  to  any  accuracy  desired  by  using 
successive  approximations. 

It  is  useful  to  observe  and  not  difficult  to  show  that  B*(a,u)  *  B*(-a,u)  and 
G(a,y)  «  -G(-a,u)  for  all  a. 

Even  though  the  function  G(a,y)  was  constructed  only  to  obtain  periodic  orbits, 
it  contains  also  information  about  the  dynamic  behavior  of  (4.1).  To  see  this,  we 
recall  that  there  is  a  two  dimensional  center  manifold  of  (4.1)  given  by  y  *  \i>(x,v). 
The  flow  on  the  center  manifold  is  given  by  the  equation 

x  «  Ax  ♦  X(x,4>(x,y),y).  (4.10) 

Any  periodic  orbit  of  (4.1)  must  lie  on  this  center  manifold  and  is  in  one-to-one 
correspondence  with  the  periodic  orbits  of  (4.10).  Formulas  (4.9)  also  give  a  one- 
to-one  correspondence  of  the  periodic  orbits  of  (4.1)  with  the  zeros  of  G(a,y).  One 
can  now  state 


Theorem  4.1.  Let  (a^.u)  satisfy  G(ag,y)  *  0  and  let  ^^(t)  be  the  corresponding 
periodic  orbit  of  (4.10).  Then  the  stability  properties  of  ^(t)  as  a  solution  of 
(4.10)  are  the  same  as  the  stability  properties  of  the  solution  aQ  of  the  scalar 
equation 

a  =  G (a, v) .  (4.11) 

Proof :  We  only  indicate  the  proof.  Since  x  «  0  is  a  solution  of  (4.10),  we  can 
introduce  polar  coordinates  x  *  (p  cos  8,-psin  8)  and  eliminate  t  to  obtain 

j|-R(e,p,y)  (4.12) 


where  R(0,O,u)  «  0,  3R(6,0,0)/3p  *  0.  One  can  apply  the  LS  method  to  this  equation 
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for  2>f-periodic  solutions  to  obtain 

,2* 

P  *  a  +  p*(a,w),  P*(a,y) (t)dt  =  0 

•*0 

1  (2v 

p*  *  R(6,a+p*,y)  -  -jjj-  I  R(6 ,a*p*  (6 j  ,y)d6 
1  f2” 

and  the  bifurcation  function  G(a,y)  =  — -  I  R(e,a+p*(a,y) (e ) ,y)d6.  The  2*-periodic 

Z”j0 

solutions  of  (4.12)  are  in  one-to-one  correspondence  with  the  solutions  of  G(a,y)=0. 
For  the  given  function  p*(a,y)  make  the  change  of  variables  p  =  r  +  p*(r,y)(9)  in 
(4.12).  Then 

ff  =  [1*1^  (r,P)(6)]'1G(r,u). 

Since  the  coefficient  of  G(r,y)  is  positive  for  r,y  small,  it  follows  that  the  sta¬ 
bility  properties  of  the  periodic  orbits  of  (4.12)  are  the  same  as  the  correspond¬ 
ing  equilibrium  point  of  r  =  G(r.y).  The  remainder  of  the  proof  follows  along  the 
same  lines  of  the  proof  of  Theorem  3.1. 

It  is  now  easy  to  give  interesting  applications  of  Theorem  4.1.  For  example, 
suppose  y  €  R  and 

Z(z,u)  =  Cj(u)t  ♦  0(|z|2) 

C(y)  =  C  ♦  Cj  (y) 


and  C(y)  has  two  eigenvalues  a(y)  ±  i8(u),  a(0)  *  0,  6(0)  =  1,  and  da(0)/dy  i  0. 
This  latter  condition  is  referred  to  as  the  Hopf  condition  and  says  that  the  curves 
(a (y )  ±  i6(u),u  small)  are  transversal  to  the  imaginary  axis  in  the  complex  plane. 
Using  the  formula  (4.7),  one  easily  shows  that 

G(0,y)  =  0,  3G(0,0)/3y  =  0,  32G(0,0)/3a3y  =  da(0)/dy. 


Thus,  the  Implicit  Function  Theorem  implies  there  is  a  unique  function  y*(a)  defined 
in  a  neighborhood  of  zero,  y*(0)  =  0,  such  that  a  1G(a,y*(a))  =  0.  Thus,  for  each 
a0>  equation  (4.1)  for  y  *  y*(aQ)  has  a  periodic  OTbit  given  by  (4.9)  with 
y  *  y*(aQ),  a  *  aQ.  The  stability  properties  of  this  orbit  are  determined  by  the 
sign  of  G(a,y*(a))  near  a^.  The  general  situation  is  the  same  as  the  one  in  Fig. 

3.1  if  one  assumes  a(y)  <  0  for  y  <  0,  a(y)  >  0  for  y  >  0.  The  result  just  stated 
is  referred  to  as  the  Hopf  Bifurcation  Theorem. 

The  Hopf  Bifurcation  Theorem  is  a  consequence  of  an  hypothesis  on  the  linear 
terms  in  z  in  (4.1).  The  specific  form  of  the  curve  y*(a)  depends  on  the  nonlinear 
terms  in  a  in  the  function  G(a,y).  Since  G(a,y)  is  odd  in  a,  the  simplest 
situation  is 

G(a,y)  *  Y0(u)a  +  Yj(u)a  ♦  ••• 

Yj (0)  i  0. 

I 

If  we  again  assume  u  is  a  scalar  and  y_(0)  =  a'(0),  where  a(y)  is  the  same  as 
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-l  2 

before,  the  curve  p*(a)  is  given  approximately  by  p*(a)  =  [a • (0) ]  Yj(0)a  ♦  ••• 
and  there  is  a  unique  periodic  orbit  near  a  -  0.  This  parabola  opens  to  the  left 
(right)  if  (o'  (0)]*^  (0)  <  0  (>  0).  Now  suppose  Re  XB  <  0  and  a' (0)  >  0.  Then 
the  solution  z  •  0  of  (4.1)  is  stable  for  p  <  0  and  unstable  for  p  >  0.  If 
Yj (0)  <  0,  then  the  periodic  orbit  exists  for  p  <  0  and  is  unstable  (this  is  called 
the  subcritical  case).  If  Yj(0)  >  0,  the  periodic  orbit  exists  for  u  >  0  and  is 
stable  (this  is  called  the  supercritical  case).  This  result  is  referred  to  as  the 
generic  Hopf  bifurcation  theorem. 

For  one  parameter  families  of  vector  fields,  the  generic  Hopf  bifurcation  is 
typical  in  the  sense  that  the  parameter  is  used  to  make  the  coefficient  of  a  in 
the  Taylor  expansion  of  G(a,p)  vanish  at  some  value  of  the  parameter.  For  two 
parameter  families  of  vector  fields,  it  is  possible  to  use  the  parameters  p  to  make 
the  first  and  third  terms  in  the  Taylor  expansion  vanish  at  say  p  =  0.  It  is  then 
natural  to  assume  that  G(a,0)  »  6a3  *  •••,  8  i  0.  The  analysis  of  the  number  of 
periodic  orbits  (two  at  most)  and  their  stability  is  easily  discussed  as  above.  The 
one  and  two  parameter  bifurcation  diagrams  for  8a3  ♦  pa  +  and  Ba5  ♦  pja3  ♦  p,+--- 
are  illustrated,  respectively,  in  Figures  4.1  and  4.2  with  the  number  of  periodic 
orbits  labeled.  The  coordinate  a  represents  the  amplititude  of  the  orbit. 

We  remark  that  the  same  type  of  proof  as  given  for  Theorem  3.1  for  the 
ordinary  differential  equation  (4.1)  holds  as  well  for  the  parabolic  equations  and 
functional  differential  equations  mentioned  in  Section  1. 


a 


Figure  4.2 

k  n  n 

5.  Homoclinic  and  heteroclinic  orbits.  Suppose  g  €  C  (RJR) ,  k  >  1,  and  the  equa¬ 
tion 

i  =  g(x)  (5.1) 

has  an  orbit  r  connecting  an  hyperbolic  equilibrium  point  xQ  to  an  hyperbolic  equi¬ 
librium  point  x  ;  that  is,  g(xQ)  *  0  =  g(Xj),  Reo(gx(xQ))  i  0,  Reo(gx(Xj))  i  0,  and 
there  is  a  nonconstant  solution  p(t)  of  (5.1)  such  that  p(t)  -»  xQ  as  t  -*>, 
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p(t)  ■*  x,  as  t  ♦  The  orbit  T  is  called  an  heteroclinic  orbit  if  xQ  i  Xj  and 

k  K  n 

an  homoclinic  orbit  if  Xg  *  Xj.  If  H  is  a  Banach  space,  h  €  C  (R  *  R  x  E,R  ), 
k  >_  1,  h(t,x,0)  *  0,  and  h  is  bounded  in  t,  then  the  problem  to  be  discussed  in 
this  section  is  to  determine  the  behavior  of  the  solutions  of  the  equation 

x  s  g(x)  ♦  h(t  ,x,y)  (5.2) 

for  t  £R  and  (x,y)  in  a  neighborhood  of  r  *  {0}. 

Since  Xq.Xj  are  hyperbolic,  there  are  solutions  Xg(t,y),  x  (t,y),  bounded 
for  t  €  R,  existing  for  y  small  such  that  xQ(t,0)  =  xQ,  Xj(t,0)  =  Xj  and  Xg(t,y), 
Xj(t,u)  are  hyperbolic.  Let  c  R  *  Rn  be  the  unstable  manifold  for  Xg(-,v) 

(Xj('.w))  and  C5^y )  the  corresponding  stable  manifolds.  For  y  =  0,  the  hypothesis 
that  T  connects  Xg  to  x^  implies  D  j^g  j*  $.  In  fact,  since  the  equation  for 
P  *  0  is  autonomous,  this  intersection  contains  an  orbit  of  (S.l).  We  want  to  deter¬ 
mine  the  nature  of  the  set  ^r.  n  S/.  for  y  i  0. 

Oy  ly 

Several  problems  motivate  this  type  of  investigation.  For  example,  if  the 
perturbation  h(t,x,y)  is  independent  of  t,  then  Xg(y),  Xj(y)  are  constant  and  one 
often  wants  to  determine  conditions  which  ensure  there  is  an  orbit  connecting  xQ(p) 
to  XjCu).  This  is  the  typical  problem  of  traveling  waves  in  parabolic  equations. 

If  Xg  »  Xj  and  the  perturbations  are  periodic  in  t  of  period  t,  and  Ut  = 

(x  :  (t,x)  €  ),  that  is  U  is  the  cross  section  of  %  at  t,  then  U  * 

uy  Xy  u,y  x>y 

Ut+T  .  Similarly,  the  cross  section  St(j  of  is  t-periodic  in  t.  If 

fl  f  ♦  ,  then  there  must  be  a  point  q  homoclinic  to  the  point  Xg(0,y)  for 
the  Poincare  map  it,  which  takes  points  in  Rn  to  the  value  of  the  solution  of  (5.2) 
at  time  t.  The  sets  u0y^S0y^  are  t*'e  sta^*^e  (unstable)  manifolds  of  Xg(0,y)  for  the 
map  it  and  have  nonempty  intersection  q.  If  they  intersect  transversal  ly  at  q,  then 
it  is  well  known  that  the  dynamics  near  q  can  be  described  by  the  left  shift 
automorphism  on  doubly  infinite  sequences  on  a  finite  number  of  symbols  (see  Section  61 

These  two  applications  are  sufficient  motivation  to  investigate  conditions  on 
h  to  ensure  that  “2^  fl  is  nonempty.  To  carry  out  this  investigation,  we  make 

a  transformation  from  the  continuous  functions  on  R  with  range  in  a  neighborhood  of 
T  to  continuous  functions  on  R  with  range  in  a  neighborhood  of  zero.  More  specifi¬ 
cally,  let 

x(t)  =  p(t+a)  +  z  (t*a) ,  a  €  R, 

and  choose  z(0)  to  be  orthogonal  to  p(0).  If  x(t)  is  a  solution  of  (5.2)  and  t  is 
replaced  by  t  -  o,  then  z(t)  satisfies 

z  *  A(t) z  ♦  f(t,z,y,a) 

A (t )  •  gx(p(t)), 

f(t,z,y,a)  •  g(p(t)  +  z)  -  gx(p(t))z  -  g(p(t))  ♦  h(t-a,p(t)  ♦  z,y) 


(5.3) 


We  consider  (5.3)  in  a  neighborhood  of  (z,v)  =  (0,0)  as  a  perturbation  of  the 
linear  equation 

(Lx) (t)  =0  Lx  =  d/dt  -  A( • ) •  (5.4) 

From  the  definition  of  A(t),  we  have  A(t)  ♦  A+  =  (Xj )  as  t  ♦  ♦»,  Aft)  •+  A  *  8X(X0) 
aj  t  +  -•  with  Reo(A+)  4  0,  Rea  (A  )  j*  0. 

Note  that  p(t)  is  a  nontrivial  solution  of  (5.4),  bounded  on  R. 

We  consider  the  operator  L  as  a  continuous  linear  operator  from  C*(F,Rn)  to 
C^flR^ln).  Our  first  objective  is  to  characterize  ./f'(L),  .  To  do  this,  the 

following  concept  is  useful. 

Let  X(t,s),  X(t,t)  =  I,  be  the  principal  matrix  solution  (solution  operator) 
of  (5.4).  Equation  (5.4)  is  said  to  have  an  exponential  dichotomy  on  an  interval  J 
with  constants  K,a  if  there  are  projections  P(s),  s  6  J,  continuous  in  s,  such  that, 
if  Q(s)  »  I  -  P(s),  then 

(i)  X(t,s)P(s)  *  P(t)X(t,s)  t,s£J. 

(ii)  |X(t,s)P(s)|  <  Ke‘a(t'S),  t  >_  s  in  J  (5.5) 

(iii)  |X(t,s)Q(s)|  <  Ke'0^'^,  s>t  in  J. 

The  operator  P(t)  is  called  the  projection  matrix  function  of  the  dichotomy. 

This  concept  is  equivalent  to  the  existence  of  a  projection  PQ  and  constants 
K,a  such  that 

|X(t,0)P0X'1(s,0)|  <  Ke*a(t‘s), 

|X(t,0)Q0X_1 ( s , 0) |  <  Ke'a(s_t), 

In  fact,  if  (5.5)  is  satisfied,  let  PQ  =  P (0)  and  observe  that  (i)  implies 

X(t,0)P0X_1(s,0)  »  X(t,0)X'1(s,0)P(s)  =  X(t,s)P(s). 

Thus,  (5.5)  implies  (5.6).  Conversely,  if  (5.6)  is  satisfied  and  one  defines 
P(t)  »  X(t,0)P0X  ^ (t , 0) ,  then  one  easily  verifies  (5.5). 

Remark .  Note  that  the  projection  P(t)  is  uniquely  determined  if  J  *  R,  but  is  not 
unique  in  other  cases.  Also,  if  J  is  finite,  then  there  is  always  an  exponential 
dichotomy  on  J.  One  can  choose  any  projection  PQ  on  Fn  and  define  P(t)  = 

Xft.OjPpX'V.O). 

Remark .  Definition  (5.5)  may  be  modified  to  apply  to  dynamical  systems  in  infinite 
dimensional  spaces  for  which  the  solution  operator  X(t,s)  is  only  defined  for  t  ^  s. 
Condition  (i)  is  replaced  by 

(i)'  X(t,s)P(s)  =  P(t)X(t,s),  t  s  in  J. 

(i)"  The  restriction  T(t,s)  <5?(Q(s)) ,  t  >_  s,  is  an  isomorphism  of 

4?(Q(s))  onto  #(Q(t))  and  we  define  T(s,t)  as  the  inverse  mapping. 


t  ^  s  in  J 
s  >  t  in  J. 


(5.6) 


The  most  interesting  cases  fot  the  interval  J  in  an  exponential  dichotomy  are 
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R+  *  [0,“),  R  ■  (-“,0]  and  R. 

Relation  (5.S)  has  a  very  simple  geometric  interpretation.  In  fact,  if  J  =  R+, 
then,  for  each  fixed  s,  there  is  a  finite  dimensional  subspace  of  Rn  given  by 
<5P(P(s))  and  called  the  stable  subspace  at  S  such  that  solutions  with  initial  values 
in  i*(P(  s))  at  s  tend  to  zero  uniformly  and  exponentially  as  t  -*  *.  If  J  =  R_ 
then  ^?(Q(s))  is  the  unstable  manifold  at  s  with  solutions  through  points  here  tend¬ 
ing  to  zero  uniformly  and  exponentially  as  t  +  The  fact  that  |P(t)|  <_  K  for 
all  t  implies  the  angle  between  the  subspaces  ^P(P(t)),  4P(Q(t))  is  bounded.  The 
angle  a(Y,Z),  0  a(Y,Z)  */2,  between  two  subspaces  Y,Z  in  Rn,  Y  fl  Z  =  {0}, 

Y,Z  +  {0}  is  defined  as 

cos  a(Y,Z)  =  sup{  :  y  £  Y  s  {0},  z  €  Z  x  (0)1. 

If  A  is  a  constant  n  *  n  matrix  with  Reo(A)  +  0,  then  x  -  Ax  =  0  has  an 
exponential  dichotomy  on  R  with  projection  P^  =  P(0)  given  by 

1  '  PA  =  271  f  C^-^dX 
■y 

where  y  is  a  closed  curve  in  C  enclosing  the  eigenvalues  of  A  with  positive  real 


We  need  several  fundamental  lemmas.  The  first  one  is  elementary  and  stated 
without  proof. 

Lemma  5.1.  Let  A(t)  be  an  n  *  n  matrix  function  defined  and  continuous  on  R.  Then 


the  equation  (5.4)  has  an  exponential  dichotomy  on  R  if  and  only  if  it  has  an 


exponential  dichotomy  on  both  (0,<»)  and  (-",0]  and  Rn  in  the  direct  sum  of  the  stable 
and  unstable  subspaces  at  zero. 

Leimaa  5.2.  Let  J  be  either  R| ,  R  or  R.  If.  *  =  A(t)x  has  an  exponential  dichotomy 
on  J  and  B(t)  is  a  continuous  n  *  n  matrix  function  on  J  with  [ B(t)  |  <5 ,  then 

x  «  (A (t)  +  B(t))x  (5.7) 

has  an  exponential  dichotomy  on  J  if  S  is  sufficiently  small. 

Sketch  of  Proof.  Consider  first  J  «  R+.  For  any  x  €  Rn  and  any  y  £  cjj(R+  ,Rn) , 
define 


(jy)  (t)  =  X(t,s)P(s)x  ♦  X(t,u)P(u)B(u)y(u)du 


j^t.i 


.00 

-  X (t  ,u 
t 


)Q(u)B(u)y(u)du. 


The  operator  is  motivated  by  the  following  consideration.  If  f  €  cjjflR^R11), 

then  one  can  show  that  the  solutions  x(t)  of  the  nonhomogeneous  equation 
x  *  A(t)x  ♦  f(t)  which  are  bounded  on  R+  must  be  given  by 

x (t )  *  X (t , s)P(s) x  ♦  [  X(t,u)P(u)f(u)du  -  f  X (t ,u)Q(u) f (u)du. 


For  5  sufficiently  small,  it  is  easy  to  show  that  the  operator  S  has  a 
unique  fixed  point  in  C^(R,Rn)  which  is  a  solution  of  (5.7)  on  R+,  is  continuous  in 
t,s,x  and,  for  each  fixed  t,s,  is  linear  in  x.  If  the  value  of  this  fixed  point  at 
t  *  s  is  denoted  by  P(s)x,  then  P(s)  :  Rn  ■*  Rn  is  a  continuous  linear  operator.  If 
X_(t,s)  is  the  principal  matrix  solution  of  (5.7),  then  the  fixed  point  of  S  is 

b 

given  by  Xg(t,s)P(s)x.  We  will  show  that  (5.7)  has  an  exponential  dichotomy  on  R+ 
with  projection  matrix  P(s). 

We  need  several  elementary  observations  to  prove  that  P(s)  is  a  projection. 
Since  P  (s)  =  P(s)  and  the  fixed  point  of  S  is  unique,  it  follows  that  P(s)P(s)  = 
P(s).  From  the  definition  of  Xg(t,s)P(s)x  and  P(s)x,  we  have 

P(s)x  =  P(s)x  -  f  X(s ,u)Q(u)B(u)Xg(u,s)P(s)x  du. 
s 

Operate  with  X(t,s)P(s)  to  obtain  X(t,s)P(s)P(s)x  =  X(t,s)P(s)x.  For  t  =  s,  this 
implies  P(s)P(s)  *  P(s)  for  all  s.  Operating  on  this  last  relation  with  P(s)  and 
using  the  fact  that  P(s)P(s)  =  £(s),  one  obtains  ¥  (s)  =  P(s)  and  P(s)  is  a  projec¬ 
tion. 

Let  Q(s)  =  I  -  P(s).  Using  the  fact  that  Xg(t,s)P(s)x  is  a  fixed  point  of  S, 
the  variation  of  constants  formula  for  Xg(t,s)  and  the  fact  that  X(t,s)P(s)  = 
P(t)X(t,s) ,  P(t)P(t)  =  P(t),P(t)[I-P(t)]  =  0,  one  obtains  P(t)XB(t,s)P(s)x  = 
P(t)Xg(t»s)x.  This  implies  P(t)Xg(t,s)Q(s)x  «  0  for  all  x.  But  then  the  fact  that 
X  (t,s)P(s)x  €  ^?(P(t))  implies  that 

b 

Xg(t,s)P(s)x  =  P(t)Xg(t,s)P(s)x  ♦  Q(t)Xg(t,s)P(s)x  = 

=  P(t)Xg(t,s)P(s)x  =  P(t)Xg(t,s)x. 

Thus,  Xg(t,s)P(s)  =  P(t)Xg(t,s)  for  all  t,s  and  (i)  in  (5.5)  is  satisfied. 

After  a  few  computations,  one  obtains 
Xg(t,s)Q(s)x 

*  X(t,s)Q(s)Q(s)x  +  |  X(t,u)P(u)B(u)Xg(u,s)Q(s)x  -  j  X(t,u)Q(u)B(u)Xg(u,s)Q(s)ds. 

If  0  =  a_1K«  <  1/2,  B  »  a(l-28)1/2,  o  =  0_1  (1- d-20)1^2] ,  then  one  can  obtain  the 
following  estimate  (it  is  nontrivial): 

|Xg(t,s)P(s)|  <  oXe"6(t's),  t  >  s  >  0. 

|Xg(t,s)Q(s)  I  <  oKe'®*5'^,  S  >  t  >  0. 

Thus,  estimates  (ii),  (iii)  hold.  This  proves  the  lemma  for  J  *  R+.  The  same  type 
of  argument  applies  to  I!_  and  R. 

Remark ■  In  the  proof  of  Lemma  5.2,  one  obtains  the  estimate 

|P(s)  -  P(s)|  <  «f,K2/(a*B) 
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which  •*  0  as  4  -»  0  uniformly  in  s. 

We  state  the  following  lemma  without  proof  since  it  is  similar  to  the  proof  of 
Lemma  5.2. 

Lemma  5.3.  If  Aft)  ■*  A~  as  t  ■»  t“,  Rea(A~)  +  0,  then  (5.4)  has  an  exponential 
dichotomy  on  R.  (R_)  with  projection  matrix  P+(t)(P"(t))  satisfying  P*(t)  +  P!  as 
t  "►  ±°°. 

Lemma  5.4.  Let  A(t)  be  an  n  *  n  matrix  function,  bounded  and  continuous  on  K  such 
that  equation  (5.4)  has  an  exponential  dichotomy  on  R|  and  R  with  projection 
matrices  P+(t).P'(t),  respectively.  Then  L  :  C*(R,Rn)  -*■  C°(R,Rn)  ^s  Fredholm  of 
index  dim  4P(P*(0))  ♦  dim  ^?( P‘(0))  -  n, 

>U)  =  {g  £  C^flRJRn)  :  g(t)-A(t)g(t)  =  0,  t£S,  g(0)  €4f(P+(0))  n^I-P'  (0))  J 

•^(L)  =  €  cj|(R,Jlr')  :  |  iji*(t)f(t)dt  =  0,  *  »  transpose. 

for  all  €  C*(*,Rn)  satisfying  the  adjoint  equation  x  #A*(t)x  =  oj. 

Remark ■  If  A(t)  satisfies  the  conditions  of  Lemma  5.3  and  PA*.PA-  are  the  PT°jectl0n 
operators  respectively  for  the  dichotomies  of  x  -  A+x  =  0,  x  -  A_x  *  0,  then 
dim^PA+)  -  dim^(P+(0)).dim  4f(PA_)  =  dim£?(P'  (0))  and  L  in  (5.4)  has  index 
dim^8tPA*)  +  dimiff PA.)  -  n. 

Sketch  of  proof  of  Lemma  5.4.  If  X(t,s)  is  the  solution  operator  of  (5.4),  then 
X  ^(t.s)  =  X(s,t)  is  the  solution  operator  of  the  adjoint  equation 

(L*x)(t)  =0,  L*  =  d/dt  ♦  A* ( • ) .  (5.8) 

This  implies  (5.8)  has  an  exponential  dichotomy  on  R  with  projection  matrix 
♦*  ♦*  -*  -* 

Q  (t)  *  I  -  P  (t)  and  on  R_  with  projection  matrix  Q  (t)  =  I  -  P  (t).  The 

fact  that  >T(L)  =  |g  €  C*flRJRn)  :  g(t)  -  A(t)g(t)  =  0,  t  £R,  g(0)  €  V  n  hj  where 

V  =  R(P+(0)),  w  *R(Q~(0)),  Q"  =  I  -  P”,  is  clear.  Also,  _^(L*)  =  |g  €  C*  (R  ,  Rn)  : 

gft)  ♦  A*(t)g(t)  =  0,  t£R,  g(0)  €  V1  n  WL|  where  V1  =  j5f(Q+*(0)),  W1  -  if(P'*(0 )). 

J  .CD 

If  f  €^f(L),  then  straightforward  calculations  show  that  i^*(t)f(t)dt  =  0 
for  all  ij/  €  _/T(L*).  Conversely,  suppose  f  satisfies  this  orthogonality  condition 
for  all  ♦  €-^tL*).  A  solution  of  Lx(t)  -  f(t),  f  £  C°(R,Rn),  is  bounded  on  R  if 
and  only  if  there  is  a  (  £Rn  such  that 

x(t)  =  X(t,0)P*(0K  ♦  joX(t,s)P*(s)f(s)ds  -  j^X(t,s)Q*(s)f(s)ds,  t  >0 

x(t)  =  X(t,0)Q'(0)£  ♦  [  X(t,s)Q'(0)f(s)ds  ♦  [  X(t ,  s)P'  (s)f  (s)ds ,  t  <0 
JO  J  -00 

that  is,  if  and  only  if 

#0  #00 

(P*(0) -Q" (0) ] C  -  X'1(s,0)P'(s)f(s)ds  *  X‘1(s.0)Q+(s)f(s)ds,  or 

>  -•  >0 


24 


o  r 

[P+(0)-Q'(0)H  =  |  p'(0)x'1(s,0)f(s)ds  +  j^Q+(0)X'1(s,0)f(s)ds. 

For  this  equation  to  have  a  solution,  one  must  have  the  right  hand  side  orthogonal 
to  all  vectors  n  €  »"  such  that  n* [P* (0) -Q' (0) ]  *  0;  that  is,  P*  (0)n  =  Q'  (0)n. 
But,  <Ji  €>f(L*)  if  and  only  if 

q-(t)  =  X*'1(t,0)(I-P+*(0))n  t  »  0 
=  X*'1(t,0)P'*(0)n  t<0 

with  P+*(0)n  *  Q**(0)n.  This  proves  that  &(L)  is  as  stated  in  the  lemma. 

Thus  L  is  Fredholm  of  index  di»>f(L)  -  dim(VL  0  KT*") .  But  dim  (V*-  f>  If*-)  * 

n-dim(V+W)  =  n-dim  V  -  dim  W  -  dim  VflK.  This  proves  the  lemma. 

We  can  now  use  Lemma  5.4  to  apply  the  LS  method  to  obtain  the  bifurcation 

equations  for  bounded  solutions  of  (5.3).  In  fact,  let  O  =  (p.d^*  ■  •  •  a  ^as^s 

for  ^(L).  Since  z(0)  is  required  to  be  orthogonal  to  p(0),  this  implies  the  pro¬ 
jection  of  t  onto  _^(L)  must  have  value  at  2ero  given  by  ®(0)b,  b  *  (0,a),  a  E 
If  we  let  t«  (♦  ,...,*)bea  basis  for  _/f(L*) ,  (L*y)(t)  *  x(t)  -  A*(t)x(t), 

ft 

f*(t)'l'(t)dt,  (I-E)f  *  CD  *  V*(t)f(t)dt  for  f€Cb(RJRn).  then  E  is  a 
od  n  q-1 

continuous  projection  on  C^CRJl")  and  4?(L)  =^(E) .  Fix  a  £R  and  apply  the  LS 

method  to  obtain  a  function  z(fi,4,e)  in  Cjj(R,Rn)  satisfying  the  equation 


z  -  A(t) z  =  Ef(-,z,u,a) 
z*(a,0.o)(0)  =  O(0)b,  b  =  (0,a). 


(S.9) 


This  function  z(a,u,a)  will  be  a  solution  of  (5.3)  if  and  only  if  (a,y,a)  satisfy 
(I-E)f(- ,z(a,u,a) ,w,a)  =  0  which  is  equivalent  to 
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(p(t),p(t))  •*  (0,0)  is  t  •*  Z«>;  that  is,  r  is  an  homoclinic  orbit  through  x  *  0. 

Since  zero  is  a  saddle  point  for  (5.12),  there  is  a  unique  hyperbolic  periodic 
solution  ip(t,u)  of  (5.11)  of  period  1  for  |w|  small,  ip(t,0)  =  0.  Our  objective  is  to 
give  necessary  and  sufficient  conditions  on  (5.11)  in  order  that  there  is  a  homoclinic 
point  to  <p(0,w)  in  a  small  neighborhood  of  (0)  *  r  for u in  a  small  neighborhood  of 
zero.  ’  From  the  definition,  we  must  determine  a  solution  x  of  (5.11)  which  remains 
in  a  small  neighborhood  of  r  for  w  near  zero  with  the  property  that  x(t)  -  <p(t,y)  ■*  0 
as  t  •*  ±«».  The  trajectory  { (t,<p(t,u)),t  £11}  being  a  hyperbolic  saddle  implies  that 
we  need  only  look  for  solutions  of  (5.11)  which  remain  in  a  small  neighborhood  of 
R  x  T  for  u  near  zero.  Thus,  for  any  solution  x(t)  of  (5.11),  we  let  x(t)  =  p(t*a)  + 
z(t+a),  replace  t  by  t  -  a  to  obtain  the  equation 


(5.13) 


22  =  -g'(p(t))z1-  u1z2-yiP(t)+'J2^(t'a)"8(P(t)  +  2j)  +  g(p(t))+g'  (p(t))Zj . 


One  can  now  apply  the  previous  theory  for  small  solutions  z(t)  of  (5.13)  for  y  in  a 
neighborhood  of  zero. 

The  linear  variational  equation  around  T  is 

Xj  *  x2,  i2  =  -g'(p(t))x2.  (5.14) 

The  only  solutions  of  this  equation  which  are  bounded  on  R  must  be  multiples  of 
(p(t).p(t)).  This  is  proved  by  using  the  fact  that  a  principal  matrix  solution  has 
determinant  equal  to  one  for  all  t.  The  adjoint  equation 

*  g'(p(t))x2,  x2  =  -Xj 

has  the  solution  (-p(t).p(t))  bounded  on  R  and  all  other  solutions  bounded  on  R  are 
multiples  of  this  one.  In  the  terminology  of  the  general  setting  for  equation  (5.1), 
(5.4),  we  have  shown  that  >fXL)  is  spanned  by  (p(t),f5(t)),  /f(L*)  is  spanned  by 
(-p(t).p(t))  and  L  has  index  zero.  Thus,  when  we  apply  the  LS  method,  q  =  p  *  1 
and  the  vector  a  is  not  needed  in  (5.9).  The  bifurcation  function  in  (5.10)  is 
given  approximately  by 


G(u,a)  *  j"°  fKtH-UjpCtj+Ujfd-cOldt  ♦  o(|y|) 
as  u  0.  This  says  that  should  be  given  approximately  by 


o  —  *  h(a) 
2 


(5.15) 


(S . 16) 


p  (t)dt,  h(a)  =  P(t)f  (t-a)dt . 


The  function  h(a)  is  periodic  of  period  1.  If  nWj/u2  is  given  and  satisfies 
minah(a)  <  nuj/ii2  f  i*«xah(a),  and  h(tjg)  *  niij /u2>  then,  if  u  is  sufficiently  small 
there  always  exist  an  a(y),a(0)  «  a^,  such  that  equation  (5.11)  has  a  solution  which 
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approaches  the  periodic  solution  <p(t,u)  as  t  ■*  ±®.  This  value  a(w)  corresponds  to 
the  initial  data  near  (p(a(M)),p(a(u)))  for  a  solution  x  of  (5.11)  which  lies  on  the 
stable  and  unstable  manifolds  of  the  periodic  solution  ip(t,w)  of  (5.11)  near  zero. 
Furthermore,  one  can  show  that  the  intersection  of  these  manifolds  is  transversal 
if  h'  (Op)  j*  0.  For  the  Poincare  map  tt,  this  implies  that  the  flow  has  a  behavior 
similar  to  the  one  shown  in  Fig.  5.1,  where  U  ,  S  are  the  unstable  and  stable  mani¬ 
folds  of  the  fixed  point  of  e  near  zero. 


In  the  next  section,  we  discuss  the  implications  of  the  existence  of  a  trans¬ 
verse  homoclinic  point. 

If  we  make  some  further  hypotheses  on  the  function  h(a)  in  (5.16),  we  can  dis¬ 
cuss  the  existence  of  homoclinic  orbits  for  a  full  neighborhood  of  u  =  (Wj.i^)  * 
(0,0).  Suppose  h(a)  has  an  absolute  maximum  (minimum)  at  <^(0^)  and 

h"^)  <  0,  h" (o^)  >  0.  (S.  17) 

Then  one  can  prove  that  there  are  (^-curves  C^,  in  u-space  with  tangents  at  u  *  0 
respectively  given  by  hfa^/n,  hfa^/n  which  divide  a  neighborhood  of  p  »  0  into 
sectors  as  in  Fig.  5.2,  such  that  there  are  no  homoclinic  orbits  in  Sj  and 

homoclinic  orbits  in  S^. 


Figure  5.2 

6.  Transverse  homoclinic  points.  To  say  more  about  the  flow  near  a  transverse 
homoclinic  point,  we  need  some  further  results  on  exponential  dichotomies,  he 
consider  only  periodic  systems 

x  *  f(t,x)  (6.1) 

where  f(t,x)  is  T-periodic  in  t.  If  it  is  the  Poincare  map  taking  points  xQ  €  Rn 
into  the  solution  through  XQ  at  time  T,  then  fixed  points  of  it  correspond  to  T- 
periodic  solutions  of  (6.1).  Let  <p(t)  be  a  T-periodic  solution  of  (6.1).  It  is 
hyperbolic  if  no  characteristic  exponents  of  the  linear  variational  equation 

x  •  fx(t,<P(t))x  (6.2) 

have  zero  real  parts.  This  is  equivalent  to  the  statement  that  no  eigenvalues  of 


sition  6.1.  If  (p(-)  is  a  hyperbolic  T-periodic  solution  of  (6.1),  then  £  is  a 


transversal  homoclinic  point  of  <p(0)  if  and  only  if  equation  (6.3)  has  an  exponen 


tial  dichotomy. 

We  also  need  the  following 

Lemma  6.2.  Fot  each  integer  V,  let  A^(t)  be  a  bounded  continuous  n  *  n  Matrix 
function  such  that  the  system 


has  an  exponential  dichotomy  on  an  interval  [t^  ^  with  constants  K,a  (indepen¬ 
dent  of  k)  and  projection  matrix  function  Pk(t).  Let  A(t)  =  Aj^t),  t  €  [t^.t^]  • 
Also,  suppose  | A^ (t) |  <,  M  for  all  t.k.  Then  there  exist  xQ  =  x0(K,a) ,  «0  *  60(IC,a) 
such  that  the  equation 

x  -  A(t)x  (6.5) 

has  an  exponential  dichotomy  on  F  if  t^  -  t^  j  >  x^  and  Ivw  -  pk-i'Wli{c 

PToof :  Only  the  ideas  will  be  given.  The  first  step  of  the  proof  is  to  construct 
n  «  n  matrix  functions  B^(t)  which  are  close  to  A^(t)  and  such  that  x  =  ^(tjx  has 
an  exponential  dichotomy  on  ^ . ty]  with  constants  3K,a  and  projection  matrix 
function  (t)  with  (tk_ j )  ■  R^.j  Ctk_j )  -  If  B(t)  -  B^ft)  and  R(t)  *  Rfc(t)  on 
[t^^j.t^],  and  Y(t,s)  is  the  principal  matrix  solution  of  x  *  B(t)x,  then 
Y(t,s)R(s)  *  R(t)Y(t,s).  One  then  shows  that  t^  -  t^_j  ^  2a’*ln3K  implies  there  is 
an  exponential  dichotomy  for  this  equation  with  constants  9K^,a/2  and  matrix  R(t). 
The  equation  (6.5)  then  can  be  considered  as  a  perturbation  of  x  =  B(t)x  to  obtain 
an  exponential  dichotomy  of  (6.5). 

Let  XjJt.s)  be  the  principal  matrix  solution  of  x  »  A^(t)x,  t^  ^  <_  t  <_  t^.  We 
construct  the  (t )  by  finding  a  nonsingular  transformation  of  variables  Sk(t)  which 
is  close  to  the  identify  such  that  Y^tt.O)  »  Sk(t)Xk(t,0)  is  a  fundamental  matrix 


of  solutions  for  the  equation  x  »  B^(t)x  on  [t^  j,tk]  an<^  let  Rk(t)  * 


^(t.OJP^COJY^  x(t,0).  The  simplest  form  for  Sk(t)  is 


,-l 


V*)  “  1  *  <Wi>  (t‘tk-I)(Jk*rI)*  Vi  1 1  i  V 

where  Jk+1  is  nonsingular.  Using  the  fact  that  (t.OJP^ (0)  *  P^ftjX^ft.O),  then 
one  sees  that  R^ft^)  *  R^jCt^)  if  and  only  if  Pv(tv-l^JV  *  JkPV-l(tk-l^  and 
Jk  is  nonsingular.  The  operator 

Jk  ■  pk(Vi>pk-i<Vi>  *  (,-Wi))(I-pk-i(Vi» 

is  a  simple  choice  for  J^.  To  show  has  an  inverse,  observe  that 

1  -  Jk  *  tvvi*  -  (l-pk(vi»HPk(tk.1)  -  Pk.j (»k. j) ] • 

Thus,  |l-Jk|  <_  1/2  if  4  *  1/4K  and  has  an  inverse. 

With  the  notation  above,  let  B(t)  «  Bk(t),  R(t)  =  Rjft),  t  €  [tk  ]»lk)  and  let 
Y(t,s),  Y(s,s)  ■  I,  be  the  principal  matrix  solution  of  x  =  B(t)x.  We  show  that 
x  «  B(t)x  has  an  exponential  dichotomy  with  constants  9K^,a/2  and  projection  matrix 
function  Q(t)  if  tkJ  -  tk  is  >_  2a  Jlog  3K.  Let  Q(0)  *  Q.  For  s  <_  t,  there  are 
integers  k  <_  j  such  that  tfc  j  <_  s  <  tk>  V  j  <_  t  <  t^  and  the  following  estimates 
hold 


,j-2. 


|V(t,s)Q(s)| 

<  |Y(t,tj_1)Q(tj.1)|'jl  |Y(ti+1.t.)(QCti)HY(tk,s)q(s)| 

<  (3K)j-k+1e'a(t-s)  9kVa(t-,)/2 


since  t  -  s  >_  t^_j  *  *k  i  (j-k-l)2a  1  log  3  K.  If  one  supplies  the  details  of  these 
computations,  the  proof  is  cooplete. 

Theorem  6.3.  Suppose  f  €  C^QR  *  Rn,Rn),  f  €  C^tR  x  i**")  and  fx  is  continuous  in 
x  uniformly  in  t,x.  For  each  integer  k,  suppose  that  the  system 

x  *  f(t,x)  (6.6) 


has  a  solution  wk(t)  defined  on  (tk.]>tk]  such  that 

x  *  fx(t,wk(t))x  (6.?) 

has  an  exponential  dichotomy  on  [tk  ^.t^]  with  constants  K,a  and  projection  matrix 
function  Pk(t)  and  the  following  conditions  are  satisfied: 

(i>  lvi(tk-ij  -  wk(Vi,l  *  5 

(ii)  IVi^-P  -pk(Vi>‘  <6 

(iii)  tk-tk^>T. 

Then  there  are  positive  constants  Cq,t0  and  a  function  6 0 ( i )  such, t ha t j_i f  t  > 


Proof:  Only  an  outline  of  the  proof  is  given.  If  a(t)  =  wk(t),  t  € 
then  w  is  continuous  except  with  small  jumps  8k  at  t^.  Lemma  6.2  implies  the 
equation  x  *  fx(t,w(t))x  has  an  exponential  dichotomy  on  R  if  t  is  large  enough 
and  6  is  small  enough.  The  next  idea  is  to  approximate  w{t)  by  a  continuous 
function  z(t)  by  linear  interpolation 

z(t)  *  w(t)  +  (tk-tk_i)_1 (t-sk)Bk_i  if  tkl  <  t  <sk 
tCO  •  w(t)  ♦  (tk-V1)-1(t-skiBk  if  sk  <  t  <  tk 

where  sfc  =  (tk_j+tk)/2.  Then  |z(t)-w(t)|  <_  6,  |  z(t)-w(t)  |  <_  t**6,  except  at  the 
points  tk-  Since  z(t)  is  close  to  w(t)  and  f^Ct.x)  is  continuous  in  x  uniformly 
with  respect  to  t,x,  the  equation  x  =  fx(t,z(t))x  has  an  exponential  dichotomy  on 
R.  One  now  considers  the  solutions  of  (6.6)  as  variations  from  z  by  letting  x(t)  = 
z(t)  +  v(t)  to  obtain 

v  «  f  (t,z(t))v  ♦  g(t,v) 

(6.8) 

g(t,v)  «  [f(t,z(t))-z(t)J  +  (f(t,z(t)+v)  -  f(t,z(t))  -  fx(t,z(t))v] . 

Now 

latt.O) I  1  |f(t,z(t))  -  f(t.w(t))|  +  |w(t)  -  z (t ) I  <  (const) 6 
except  at  the  points  tk<  Also, 

Igyit.v)!  *  |fx(t,z(t)+v)  -  fx(t,z(t))|  <_  U  (  |  V I  ) 

where  w  is  the  uniform  modulus  of  continuity  of  f  .  he  have  assumed  that  u(s)  ■*  0 
as  s  -*  0.  Since  v  =  fx(t,z(t))v  has  an  exponential  dichotomy  on  R,  the  equation 

v  =>  f x (t,z(t))v  +  g(t) 

where  g  is  bounded  on  R  and  continuous  except  at  the  tk>  has  a  unique  solution  Kg 
bounded  on  R  and  sup|(Kg)(t)|  <_  (const)  sup|g(t)|.  If  we  let  F  :  R  *  C°(R,Rn)  •* 

0  -  t  t  b 

)  be  defined  by  the  relation  F(t,v)  =  v(t)  -  K(g(  •  ,v(  • ))  (t) ,  then  equation 

(6.8)  has  a  solution  in  C^(RJln)  if  and  only  if  F(-,v)  =  0.  The  function  F(-,v)  is 
continuous  together  with  its  first  derivative  DvF(,,v)  in  a  neighborhood  U  of  v  =  0 
since  fx(t,x)  is  continuous  in  x  uniformly  in  (t,x).  Furthermore,  DyF(-,0)  =  I. 
Also,  F(-,<p)  «  0,  where  <p  *  K{f(',z('))  -  *(•)]•  Since  |ip|  <_  (constant)  -5,  the 
Implicit  Function  Theorem  will  imply  that,  for  6  sufficiently  small,  there  is  a 
unique  solution  v*  €  C^0RJRn)  of  F(-,v)  =  0  in  a  neighborhood  of  v  =  0  which  is 
0(|<p|)  as  |tp|  *  0.  The  function  x*(t)  *  v*(t)  ♦  z(t)  is  continuous  and  satisfies 
x*(t)  *  f(t,x*(t))  except  perhaps  at  t  >  t,.  But,  since  f(t,xl  is  continuous,  it 


and,  in  addition,  is  T-periodic  in  t.  Suppose  there  is  a  doubly  infinite  sequence 
(u^(t)}  of  hyperbolic  T-periodic  solutions  of  (6.6)  and  another  sequence  <vk(t)} 
of  bounded  solutions  such  that  vk(t)  -  u^  ^(t)  ->0  is  t  «  v^(t)  -  u^Ct)  -►  0 

as  t  ■*  and  the  equation 

x  *  fx(t,vk(t))x  (6.9) 

has  an  exponential  dichotomy  on  R  with  constants  K,a  incependent  of  k.  Then  there 
are  >  0  and  a  function  Mg(e)  such  that,  for  any  0  <  e  <  and  any  positive 
integer  m  >  M^s),  system  (6.6)  has  a  unique  solution  x(t)  defined  on  R  satisfying 

| x(t* (2k-l)mT)  -  vk(t)|  <  c 

for  -mT  5.  t  £  mT  and  all  k. 

Proof.  Only  an  outline  of  the  proof  is  given.  Let  Pj^ ft )  be  the  projection  matrix 
function  for  the  dichotomy  of  (6.9).  Since  vk(t)  -  u^t)  ♦  0  as  t  *  »,  the 
equation  *(t)  *  f^t.UjJOJx  has  an  exponential  dichotomy  on  [0,»).  Since  ^(t)  is 
periodic  in  t,  the  Floquet  theory  implies  this  equation  has  an  exponential  dichotomy 
on  R.  Let  Q^ft)  be  the  corresponding  projection  matrix  function.  An  extension  of 
Lemma  5.3  implies  I  (t) -Qk  (t)  |  +  0  as  t  +  •>.  |?k(t)-Qk_j  (tj  |  -  0  as  t  •*  - 
uniformly  with  respect  to  k.  Let  tfc  «  2kmT,  wfc(t)  «  vk(t-(2k-l)mT),  Pk ft)  = 

P(t- (2k-l)mT) ,  and  apply  Theorem  6.3  to  complete  the  proof. 

Figure  6.1  should  assist  the  reader  in  understanding  the  meaning  of  the 
hypotheses  in  Theorem  6.4  and  also  to  feel  intuitively  why  the  conclusion  is  true; 
that  is,  how  one  should  be  able  to  switch  from  comparing  x(t+(2k-l)mT)  to  vk(t)  on 
[-mT.mT]  to  comparing  x(t*(2(k+l)-l)mT)  to  vk+J(t)  on  [-mT.mT];  that  is, 
x(t+ (2  k-l)mT)  to  vk+j(-2mT*t)  on  [-mT.mT],  The  hypotheses  imply  that 
vk+l  (*2mT+T)  |  t>bt  *  vjc+j(**'n  is  c^ose  t0  CmT)  if  m  is  sufficiently  large.  Thus, 
if  x(t+ (2k-l)mT)  is  close  to  vk(t)  on  [-mT.mT],  then  x(mT+ (2k-l)mT)  = 
x(-mT+(2(k*l)-l)raT)  is  close  to  vk+1(-mT). 


Figure  6.1 


We  now  use  Theorem  6.4  to  obtain  information  abi 


tho  flow 


near  a 


t  r\o:  ?■  ver?t- 


MP 


homoclinic  point.  let  N  he  a  positive  integer  and  let  s^  Ne  tee  -et  •*'  .ili  t-.n.. 

infinite  sequence*  a  =  ; a  ,,a  ,  ,:i  ,a,,a . ■  nth  a  £  •  ■  • ,  1 . Sir  u<  j 

put  the  product  topology  on  S^,,  that  is.  a  neighborhood  has:.-  of  a  point  t 
(.  .  .  ,b  ,,b  ]  •  >  b  j  it, ,  .  .  .  1  consists  of  the-  sets  IJi  =  la  £  S  a  -  *"  for  V 

The  (right)  Bernoulli  shift  ?  of  is  defined  as  (ria  ■■  »  3  ,  . 

Corol  iary  6.5.  f  Shad  owing _ lemma)  .  Suppose  f't,x>  satisfies  tjhc  snoot  h;,e  ■><  proper;  ics 

in  Theorem  6.  3,  is  T -pcrio.il  c  in  t .  has  an  hyperbolic  1  -periodic _ <ol-.it  t„r.  u;t  1  and 

another  solution  v (. t )  such  that  Vfti  -  u(t)  -»  (I  as  t  -  t»  and  the  equation. 

x  »  f  !  t ,  v  ( t  i )  x 
x 

has  an  exponential  dichotomy  on  H  (that  is,  vIO  is  transverse  hi.moclinic  to  u!0!  bo 
Proposition  6 . 1 ) . 

Then  there  exists  an  t  >  0  and,  for  each  positive  integer  S,  a  fun. tier 
M^.fe)  such  that,  for  any  o  •  t  <_  t(1  and  any  positive  integer  r  _  V  i-  I  and  ary 
a  £  S^,  equation  (o.o'i  has  a  unique  solution  x  i t )  satijsfVing 

jx  ft*!  Jk  - 1  )  m  l'  1  -  v  f  t  *a,  'I  i  i  - 
a  k 

for  -mT  ^  t  <_  mT  and  all  k . 

The  mapping  q>!a)  -  x^fo)  is  a  homeomorphi sn  ont  v  a  compact  set  of  lRn  cn  which 
the  iterate  r'"1  of  the  Poincare  map  -  is  invariant  and  -~w  op  =  P  of.,  where  £ 

is  the  right  Bernoulli  shift  on  S  . 

Proof :  Let  v^(t)  =  v(t*a^T),  u^Ct)  =  uft),  k  -  0,  1....  and  apply  Theorem  6.4  to 
obtain  the  existence  of  x  ft).  Using  the  uniqueness  of  x  (t),  one  can  prove  after 
some  computations  that  ipfal  is  a  homeomorphi  sm.  To  show  that  tt“  o  P  =  P  o  S,  ob¬ 
serve  that 

|xa(t*(2(K*l)-l)mT)  -  v(fak+1T)|  *  e 

for  -mT  <  t  <  mT  for  all  t.  Thus,  uniqueness  implies  xa(t*2mT)  =  Thus, 

P(8(a))  =  xa(2mT)  =  *2m(V>(a)). 

Remark .  The  "shadowing  lemma"  for  dif feomorphisms  is  usually  proved  by  using  horse¬ 
shoes.  In  this  type  of  proof,  the  symbols  (0,1 . N-l)  occur  as  a  specification  of 

whether  or  not  iterates  of  a  point  belong  to  certain  intervals. 

7.  Codimension  one  bifurcations  in  the  plane.  In  previous  sections,  we  have  dis¬ 
cussed  various  types  of  dynamic  bifurcation  for  autonomous  systems,  for  example,  the 
saddle-node  bifurcation  in  Section  3  and  the  generic  Hopf  bifurcation  in  Section  4. 

Ne  also  discussed  some  aspects  of  homoclinic  bifurcation,  lor  differential  equations 


a. 
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in  the  plane,  much  more  information  is  available.  In  fact,  one  can  completely 
characterize  all  of  the  codimension  one  bifurcations.  To  make  this  more  precise, 
we  need  the  concept  of  structural  stability.  We  restrict  the  discussion  to  the 
interior  S  of  a  closed  curve  r  without  contact  to  any  of  the  vector  fields  to  be 
considered.  Let  the  set  of  all  such  Cr  vector  fields.  Two  vector  fields 

X,Y  in  r  >_  1,  are  equivalent  if  there  is  a  homeomorphism  on  SI  U  r  which  maps 

orbits  of  one  onto  orbits  of  the  other  and  preserves  the  sense  of  direction  in  time. 
This  is  an  equivalence  relation"’-'1  among  vector  fields.  X  is  structurally  stable  if 
every  Y  in  a  neighborhood  of  X  is  equivalent  to  X. 

The  condition  that  the  vector  fields  are  nowhere  tangent  to  r  is  very  con¬ 
venient  since  it  makes  the  domain  where  the  differential  equations  are  being  con¬ 
sidered  to  be  a  compact  set.  It  also  avoids  certain  complications  which  can  arise 
at  the  boundary.  Our  hypothesis  does  put  restrictions  on  the  vector  fields.  Since 
our  objective  is  to  present  some  of  the  basic  ideas,  the  hypothesis  seems  justified 
in  view  of  the  technicalities  that  arise  in  considering  the  noncompact  case  or  flows 
on  manifolds. 

The  basic  result  on  structural  stability  is  the  following. 

Theorem  7.1.  An  f  is  structurally  stable  if  and  only  if  every  equilibrium 

point  and  every  periodic  orbit  is  hyperbolic  and  there  are  no  connections  between 
saddle  points.  The  set  of  structurally  stable  systems  is  obviously  open  but  is  also 
dense  in  <^j. 

An  X  £2^  is  a  bifurcation  point  (a  vector  field  for  which  a  perturbation  could 
lead  to  a  bifurcation)  if  X  is  structurally  unstable;  that  is,  not  structurally 
stable.  We  now  give  an  inductive  definition  of  a  bifurcation  point  of  codimension  k. 

The  vector  field  X  is  a  bifurcation  point  of  codimension  0  if  it  is  structurally 
stable.  X  is  a  bifurcation  point  of  codimension  1  if  it  is  not  codimension  zero  and 
there  is  a  neighborhood  of  X  which  has  only  bifurcation  points  of  codimension  0  or 
ones  which  are  equivalent  to  X.  It  is  a  bifurcation  point  of  codimension  2  if  it 
is  not  of  codimension  zero  or  one  and  there  is  a  neighborhood  containing  only  bifur¬ 
cation  points  of  codimension  0  or  1 ,  or,  ones  which  are  equivalent  to  X.  Similarly, 
one  defines  bifurcation  points  of  codimension  k. 

The  following  result  is  a  classification  of  bifurcation  points  of  codimension 
one  in  the  plane. 

Theorem  7.2.  A  vecotr  field  f  €2^.  r  ^  3,  is  a  bifurcation  point  of  codimension  1 
if  and  only  if  theTe  is  a  neighborhood  W  of  f  and  a  submanifold  T  of  codimension 
one  in  W  such  that  W  r  *  Uj  U  where  each  g  €  U.  is  structurally  stable  but 

g  4-  h  ijf  g  €  Uj,  h  €  .  For  g  €  r,  only  one  of  the  following  situations  prevails: 

(i)  g  €  T  has  an  elementary  saddle-node  at  x(),  there  are  no  euui  librium 


points  of  g  near  x^  ijf  g  £  Uj  and  a  saddle  and  node  near  x^  j_f  g  €  U^. 

(ii)  g  €  T  has  an  elementary  focus  at  xQ,  there  is  no  periodic  orbit  of  g 
near  x^  i_f  g  £  Uj  and  a  periodic  orbit  near  if  g  €  --the  generic 
Hopf  bifurcation. 

(iii)  g  €  r  has  a  periodic  orbit  y  which  is  stable  from  one  side,  unstable 
from  the  other,  g  €  Uj  has  no  periodic  orbit  neaT  y  and  g  £  has 
two  hyperbolic  periodic  orbits  near  y. 

(iv)  Op  =  tr  3f(x0)/3x  +  0,  g  €  r  has  a  homoclinic  orbit  containing  a 

saddle  point  x^,  g  €  has  a  saddle  near  xQ  and  no  periodic  orbit 
near  y,  g  £  has  a  saddle  point  and  a  unique  hyperbolic  periodic 
orbit  near  y  which  coalesce  as  g  r. 

(v)  there  is  a  connection  between  distinct  saddle  points. 

Each  of  the  cases  (i)-(v)  is  shown  in  Figure  7.1. 


u,  r  u2 


Figure  7.1 

We  do  not  give  a  proof  of  either  Theorem  7.1  or  7.2.  We  only  remark  that  we 
have  discussed  the  bifurcations  that  occur  in  cases  (i)  and  (ii)  in  Section  3  and 
4.  Case  (iii)  can  be  discussed  using  the  methods  in  Section  4  after  the  introduc¬ 
tion  of  a  coordinate  system  x  t— *  (6,p),  x  *  p(o)  ♦  pv(«)  where  p(t)  is  a  non- 
hyperbolic  periodic  solution  of  a  bifurcation  point  f  e  and  v(o)  is  orthogonal 
to  p(o)-  Case  (iv)  is  the  most  difficult.  Let  us  only  motivate  the  condition  n  /  0. 
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Let  us  suppose  that  we  have  a  vector  field  f(x,p)  in  R*"  depending  on  a  scalar 
paraaeter  p  with  the  property  that  there  is  a  hyperbolic  saddle  point  at  0  with  the 
eigenvalues  bounded  away  from  zero  for  all  p.  Also,  suppose  there  is  a  hyperbolic 
periodic  orbit  Py  of  period  w(y)  which  has  the  property  that  dist(0,Py)  0  as 
p  ■*  0.  This  creates  a  homoclinic  orbit  T  at  p  =  0.  For  such  a  situation  to  occur, 
the  period  w(u)  must  approach  “  as  p  •+  0.  If  this  one  parameter  family  of  vector 
fields  is  generic,  then  one  cannot  expect  to  have  other  bad  things  happen  as  we 
change  the  vector  field  since  we  use  the  parameter  to  make  d(p)  -*■  «  as  p  ■*  0.  In 
particular,  the  rate  of  attraction  or  repulsion  of  the  periodic  orbit  Py 
should  be  exponential  and  uniform  in  p.  If  we  keep  this  uniformity  as  p  +  0,  then 
the  orbit  T  should  have  the  property  that  it  is  either  asymptotically  stable  or 
unstable  depending  on  whether  Py  is  stable  or  unstable. 

To  find  a  quantitative  expression  for  this  uniformity,  let  us  recall  the 
formula  for  the  characteristic  exponents  of  the  linear  variational  equation  for  a 
periodic  orbit  in  the  plane.  If  P  ■  {p  (t)  :  0  t  <  w(v)}  and  x  *  py  ♦  y,  then 
the  linear  variational  equation  for  py  is 

y  =  5;  cpp<t>.")y- 

The  nontrivial  periodic  function  py  satisfies  this  equation  and  thus  one  charac¬ 
teristic  exponent  may  be  taken  to  be  zero.  Let  X(p)  be  the  other  characteristic 
exponent.  The  sum  X(p)  ♦  o  =  X(u)  of  the  characteristic  exponents  must  be 

,  f“(R) 

=  J0  (tr  9f(PpCt).P)/dx]dt. 

Then  one  can  show  that 

,  f“(R)/2 

*oo  =  ^nrr  ttr  9f(pu(t).p)/9x]dt-  tr  af(o,o)/ax 

1  J  1  -iu(u)/2  M 

as  w  -*■  0.  Consequently,  the  rate  of  attraction  or  repulsion  of  each  Py  will  be 
exponential  and  uniform  in p if 

def  afto^oi  =  0> 

0  3x 

The  fact  that  only  two  possibilities  arise  in  a  neighborhood  of  a  bifurcation 
point  of  codimension  one  suggests  that  this  is  the  typical  or  generic  situation  that 
arises  in  the  discussion  of  one  parameter  families  of  vector  fields.  This  is,  in 
fact,  the  case  and  one  can  prove 

Theorem  7.3.  if  ®k  =  f«P  :  [0,1]  *3^,  q>  €  Ck },  k  >  3,  and  =  {<P  €  :  q>(t)  is 

structurally  stable  except  at  a  finite  number  (depending  on  <P)  points  t  with  «p(t.) 

*  k  ■  k  * 

a  bifurcation  point  of  codimension  one) ,  then  ®  is  a  residual  set  in  «  . 


The  analysis  of  bifurcations  of  codimension  greater  than  one  are  generally 
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auch  more  difficult.  For  some  cases,  the  ideas  are  clearly  understood  and  it  is 
mainly  a  technical  problem  to  do  the  complete  analysis.  This  remark  applies  to  the 
situation  where  the  linear  variational  equation  near  an  equilibrium  point  has  either 
a  simple  eigenvalue  zero  or  a  pair  of  purely  imaginary  eigenvalues  on  the  imaginary 
axis.  The  ideas  in  Sections  3  and  4  apply  to  this  situation.  In  the  plane,  the 
methods  necessary  to  analyze  the  bifurcations  of  higher  codimension  resulting  from 
the  nonhyperbolicity  of  a  periodic  orbit  are  also  clear.  For  other  situations, 
special  difficulties  arise  and  each  problem  is  a  challenge  in  itself.  In  the  next 
sections,  we  discuss  some  special  codimension  two  bifurcations. 

8.  Two  zero  eigenvalues.  In  this  section,  we  discuss  a  codimension  two  bifurcation. 
The  unperturbed  system  is  taken  to  be 

x  *  -y,  y  =  ax2  +  6xy  (8.1) 


where  a*  0,  g  i  0.  Without  loss  in  generality,  one  can  assume  a  <  0,  6  >  0.  The 
perturbed  system  will  be 

•  •  2 

x  «  y,  y  =  EjX  ♦  c^y  ♦  ax  ♦  6xy  (8.2) 

where  Cj.e^  are  SBa11  parameters.  The  problem  is  to  discuss  the  behavior  of  the 
solutions  of  (8.2)  in  a  neighborhood  of  (x,y)  *  (0,0)  for  (Cj.e^)  in  a  neighborhood 
of  (0,0).  We  remark  that  the  conclusions  below  are  valid  for  some  higher  order  per¬ 
turbations  of  (8.2)  and  that  we  omit  such  terms  only  for  simplicity  in  notation. 

We  only  discuss  Cj  >  0,  since  the  other  case  is  less  intersting.  The  scaling 
2  2 

tj  *  4  ,  e2  *  u*  >  6  >  0 

t. — ►«'1t,  x  i— *  «2|a|’1x.  y  w—w  63|a|'1y 


leads  to  new  equations 

.  2 

x*y,  y-x-x  ♦  uiy  *  4yxy  (8.3) 

where  y  »  B|a|  . 

For  4*0,  Equation  (8.3)  becomes  the  conservative  system 


with  first  integral 


x  *  y,  y  *  x 


x 


2 


V(x,y) 


(8.4) 


(8.5) 


The  equilibrium  point  (0,0)  is  a  saddle  with  a  homoclinic  orbit  through  it  while  the 
equilibrium  point  (1,0)  is  a  center. 

The  parametri zation  of  the  phase  space  by  the  above  scaling  is  suggested  by 


the  following.  The  orbits  of  the  unperturbed  equation  (8.1]  satisfy  the  equation 
ydy  *  (ax  +gxy)dx.  The  orbit  which  passes  through  the  origin  is  given  approximately 
by  y  =  ax  .  Thus  when  one  "blows  up"  the  flow  at  the  origin,  it  is  natural  to  do 
it  by  parametrizing  the  phase  space  with  cusps.  The  orbits  for  (8.2]  satisfy 

dy  *  I(£j  ♦  ax]xy_1  +  (E2+8x)]dx. 

For  >  0,  it  turns  out  that  the  appropriate  parameterization  in  parameter  space  is 
to  take  Cj.e^xofthe  same  order.  This  leads  to  the  scaling  used  above.  Notice  that 
the  flows  defined  by  (8.2]  and  (8.3]  are  equivalent  for  6*0,  but  are  not  equivalent 
for  6=0. 

Our  next  task  is  to  determine  the  curves  in  the  (Cj.O-plane  (that  is,  the 
values  of  6,p)  at  which  the  topological  structure  of  the  traiectories  of  Equation 
(8.3)  changes.  As  we  shall  see,  this  structure  can  change  only  due  to  a  homoclinic 
orbit  or  a  change  in  stability  of  the  equilibrium  point  (1,0). 

Let  r  *  ((q(t),q(t)),  t  €  ®]  U  {(0,0)}  be  the  homoclinic  orbit  of  (8.4).  he 
may  now  apply  the  results  of  Section  S  to  obtain  the  curve  in  (u, 6) -space  such  that 
(8.3)  has  a  homoclinic  orbit.  In  particular,  from  Eq.  (S.4),  we  have  this  curve 
corresponds  to  the  zeros  of  the  bifurcation  function 

G(u,6)  *  y  ♦  yv  ♦  G (y,6) 


r  •  2 ,  r 

*  MW/ 

J  -  oo  J  -a 


where  G(p,0)  =  0.  There  is  no  a  in  G(p,6)  since  the  equation  is  autonomous.  One 
can  show  that  v  *  6/7.  The  equation  G(p,6)  =  0  has  a  unique  solution  (u(6),6), 

0  i.  I4I  i  so’  40  *  *  u0  =  ’YV‘  Finall>'’  the  curve  in  (ej , e^) -space  along 

which  there  is  a  homoclinic  orbit  is  given  by 

C«,  *  {(tj.Cj)  :  C2  *  v(Cj^^)  E  j  ,  v(0)  =  -yv)  . 

On  the  curve  x  {0,0}  we  have  °q  *  e2  <  0,  where  Og  is  the  number  given  in  Theorem 
7.2.  From  part  (iv)  of  Theorem  7.2,  this  suggests  there  should  be  a  periodic  orbit 
near  this  curve. 

We  now  discuss  the  periodic  orbits  of  (8.3).  This  part  of  the  analysis  in 
problems  of  this  type  is  the  most  difficult,  especially  the  discussion  of  the  number 
of  periodic  orbits  that  can  exist.  One  can  prove  the  following  lemma. 

Lemma  8 . 1 .  Every  periodic  orbit  of  Equation  (8.3)  must  intersect  the  segment 
(0,1)  x  (0)  in  the  (x.y)-plane.  There  is  a  continuous  positive  function  : 


f0-1)  -»*  and  a  continuously  differentiable  function  u*(b,6),  b  e  (0,1),  | « |  <  *  (h) 
such  that  there  is  a  periodic  orbit  of  Equation  (8.3)  through  (b.O)  if  and  only  if 
v  *  u*(b,6).  Furthermore, 
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a(b) 


u*( b.O)  =  -YS(b)/a(b) 

c(b)  fC  (b) 

ydx,  6(b)  =  xydx  >  0 

b  'b 


y  >  o, 


id 

2 


2  K2 

x  -b 


3  .3 
x  -b 


c(b)  >  1, 


2  2  3  3 

c  (b)  V  c  (b)-b 

2  '  3 


and  du*(b,0)/db  <  0.  Also,  u*(b,0)  +  -y  as  b  +  1,  w*(b,0)  -*•  -  yy  as  b  +  0.  Finally, 
if  v  -  u*(b,5)  for  a  fixed  b  6  (0,1)  and  j <5 1  <  ^(b),  then  the  periodic  orbit  through 
(b.O)  is  the  only  one  corresponding  to  this  y,6. 


Remark :  The  assertion  dw*(b,0)/db  <  0  implies,  for  any  bQ,  there  is  a  6Q (bQ)  such 
that  along  the  curve  Cj  =  u*(b0>6)e2,  0  <  6  <  6Q (b) ,  there  is  a  unique  periodic  orbit 
of  (8.3)  which  approaches  the  periodic  orbit  of  (8.4)  through  (bQ,0)  as  5  -*•  0. 

We  give  only  an  indication  of  the  proof  of  Lemma  8.1.  If  V(x,y)  is  the  deriva- 

2 

tive  of  V  along  thi  solutions  of  (8.3).,  then  V  «  S(u+yx)y  . 

For  a  fixed  (b.O),  0  <  b  <  1,  and  6  sufficiently  small,  there  are  numbers 
tj  =  Tj(b,6,u)  <  0  <  t2  *  t2(b,6,u)  such  that  the  solution  through  (b.O)  intersects 
the  x-axis  at  time  tj  at  a  point  larger  than  1.  Furthermore,  for  t  €  (ty,^),  the 
orbit  inter?ects  the  x-axis  only  for  t  =  0.  Let  r  *  T (b , 6 , u )  be  that  part  of  the 
orbit  through  (b,0)  corresponding  to  t  €  [tj.tj].  For  r  to  be  a  periodic  orbit,  it 
is  necessary  and  sufficient  that  Vdt  «  0,  which  for  5  i  0,  is  equivalent  to 

Jr 

F(b,«,u)  d£f  [  (u+yx)y2  =  0. 

h 

One  can  now  apply  the  Implicit  Function  Theorem  to  this  equation  near  the 
point  (bg.O.Ug),  *  -y8(b0)/a(bg) .  This  will  prove  the  first  part  of  the  lemma. 

The  fact  that  p*(b,0)  approaches  the  limits  indicated  above  require  only  elementary 
computations. 

The  difficult  part  of  the  lemma  is  to  show  du*(b,0)/db  <  0.  If  2a  =  -b"  * 
2b3/3,  then  0  <  b  <  1  implies  -1  <  6a  <  0.  Let  b(a)  be  the  inverse  of  this  trans¬ 
formation  and  put  v(a)  =  6(b(a)/a(b(a)) .  If  we  consider  6,  as  functions  of  a 
and  let  "  '  "  be  differentiation  with  respect  to  a,  then  the  lemma  is  proved  if 
one  shows  that  v'  <  0,  -1  <  6a  <  0.  To  carry  out  this  proof,  one  exploits  special 
properties  of  the  elliptic  integrals  a, 6.  More  precisely,  one  shows  that  a, 6  and 
a", 6"  can  be  expressed  as  linear  combinations  of  a’, 6'.  We  are  going  to  show  that 


Sa  *  6aa'  *8'  ,  358  e  6aa'  ♦  6(l*5a)8' 

6a(l+6a)a"  *  -6aa '  -  6'  ,  (l+6a)8"  =  6'  -  a'  . 


(8.6) 
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We  have  already  proved  that  8' (a)  =  j  xy^dx.  Using  the  fact  that  x  * 

yy  ♦  x2,  y(b)  =  y(c)  *  0,  one  obtains  6'  (a*?  =  x2y*3dx.  Integrating  by  parts, 

x  2  J  b 

using  yy^  *  x  -  x  and  the  formula  for  8‘,  we  have 

a  =  [  ydx  =  - f  xy  dx  *  - (  x(x-x2)y_1dx  =  -6'  +  f  x3y_1dx.  (8.7) 

Jb  '  b  x  Jb  Jb 

This  relation  and  the  formula  for  y  imply  that 

fC  .  fc  2  -1  .  fc,  2  2  3  ,  ,  -1. 

a  =  ydx  =  I  y  y  dx  =  [x  -  -rx  ♦  2a]y  dx 

J b  Jb  J  b  3 

*  8’  -  |(a+8')  +  2aa* 


or, 


5a  =  6aa '  +  8 ' . 


Integrating  8  by  parts  and  using  the  formula  for  yy^,  we  have 


Txydx  »  -p 

Jh  Jh 


(xV2)yxdx  = 


lfc  3  -1.  lfc  4  -l 
jj«,  dx.jjx,  dx. 


Using  (8.7),  this  implies  8  ♦  (1/2)  (a+S* )  *  (1/2) 
2 

y  ,  (8.7)  and  this  latter  relation,  one  obtains 


fC  4  -1,. 

x  y  dx. 
Jh 


Using  the  formula  for 


8*|  xy 2y  *dx  =  j  (x3-  |x4  ♦  2ax]y‘1dx  =  a  ♦  6'  -  -  |(a+6')  ♦  2a6'  . 

Simplifying  this  expression,  one  obtains  76  *  a  +  8'  ♦  6aB' .  Using  the  previously- 
obtained  expression  for  a  in  terms  of  a', 8',  one  finds  that  358  *  6aa’  ♦  6(l+5a)8'. 
The  expressions  for  a", 8"  are  obtained  from  the  relations  for  a, 8.  This  completes 
the  proof  of  (8.6). 

Using  these  relations,  one  now  proves  that,  if  v'(a)  *  0  for  some  a, 

-1  <  6a  <  0,  then 

-6a(l+6a)av"/a'  «  -  (v  ♦  6a)2  -  6a(l*6a)  <  0; 


that  is,  v"(a)  <  0. 

Next,  one  shows  that  v'  (a)  *  0  implies  7v2(a)  +  6(2a-l)v(a)  -  6a  =  0.  This 
implies  that,  if  v'(a)  *  0  and  v(a)  *  1,  then  6a  =  -1.  Since  v(0)  =  6/7,  we  have 
v' (a)  *  0  implies  0  <  v(a)  <  1. 

Using  the  fact  that  Uq(-1/6)  *  1,  Wq(0)  *  6/7,  one  easily  concludes  that 
y'(a)  <  0,  -1  <  6a  <  0  and  the  lemma  is  proved. 

Using  Lemma  8.1  and  the  remarks  about  the  homoclinic  orbit  before  the  state¬ 
ment  of  the  lemma,  Theorem  7.2  (iv)  implies  that  for  each  point  in  the  region  below 
the  curve  C^,  there  is  a  unique  periodic  orbit.  Next,  we  analyze  the  behavior  of 
the  solutions  of  (8.3)  near  the  equilibrium  point  (1,0).  This  point  is  a  stable 
focus  if  u  <  -y,  and  an  unstable  focus  if  u  >  -y,y  *  6 | a [  *.  The  curve  u  *  -y  is 
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therefore  a  possible  value  for  a  Hopf  bifurcation.  One  can  apply  the  method  of 
Section  4  for  the  periodic  orbits  near  (1,0)  and  obtain  a  bifurcation  function 
G(a,u,6)  for  | £ |  <  6g,  |  aj  <  aQ,  |v  +  y|  <  n  for  some  constants  6Q  >0,  aQ  >  o, 
n  >  0.  Since  Equation  (8.3)  for  6=0  has  a  center  at  (1,0),  it  follows  that 
G(a,u,0)  =  0  for  all  a,ji.  Thus,  the  appropriate  bifurcation  function  to  consider 
is  H(a,u,6)  *  G(a,y,6)/6.  This  function  satisfies  H(0,w,6)  =  0,  3H(0,(j ,0)/S a  = 

(u  ♦  r)/2.  Lemma  8.1  implies  for  each  u  >  -y  and  sufficiently  close  to  -y,  there 
is  a  unique  periodic  orbit  of  (8.3)  through  (b,0)  with  b  near  1.  Thus,  if  b  is 
taken  close  enough  to  1,  this  periodic  orbit  must  correspond  to  a  zero  of  the 
bifurcation  function  H(a,y,6).  This  proves  there  is  a  Hopf  bifurcation  at  u  =  -y 
and  the  periodic  orbit  is  asymptotically  stable. 

One  thus  obtains  the  complete  bifurcation  diagrams  as  shown  in  Fig.  8.1  with 
the  flow  in  each  sector  given  in  Fig.  8.2.  We  draw  the  curves  in  Fig.  8.1  as 
straight  lines  but  this  is  really  only  the  first  approximation. 


Figure  8.2 

9.  Two  zero  roots  with  symmetry.  In  a  planar  system  where  the  matrix  of  the  linear 
variational  equation  near  an  equilibrium  point  has  both  eigenvalues  zero  with  non- 
simple  elementary  divisors,  the  analysis  in  the  previous  section  showed  that  the 
vector  field  (8.1)  with  quadratic  terms  was  a  codimension  two  bifurcation.  If  there 
is  some  symmetry  in  the  vector  field;  for  example,  it  is  odd  in  (x,y),  then  the 
quadratic  terms  in  the  Taylor  expansion  vanish.  Thus,  it  becomes  of  interest  to 
know  what  additional  nonlinear  terms  are  needed  in  order  to  obtain  a  codimension  two 
bifurcation.  In  this  section,  we  summarize  some  results  with  only  brief  indications 
of  the  proofs. 


Consider  the  equation 

3  a  2 

x  «  y,  y  =  CjX  ♦  e2y  ♦  ax  ♦  6x  y 


with  a  <  0,  B  <  0  and  tj.c^  *mal*  parameters.  The  problem  is  to  analyze  the 


(9.1) 


behavior  of  the  solutions  of  (9.1)  in  a  neighborhood  of  (x,y)  =  (0,0)  for  (cj,e2) 
in  a  neighborhood  of  (ej>e2)  1  (0,0). 

The  bifurcation  diagram  is  shown  in  Fig.  9.1  with  the  flow  in  each  sector 
given  in  Fig.  9.2. 


REGION  I  REGION  Z 


Figure  9. 


Let  us  give  an  idea  of  how  these  results  are  obtained  for  the  case  Ej  >  0. 
Firstly,  one  introduces  scaling 

*  =  I £ja_1 1 1/2.  e2  =  I “ 1 1/2*tj »  *  1 — *  6x»  y  «2M1/2y 

t l.r^V1  t 


to  obtain 


where  y  *  sla 


•  •  3  2 

x  =  y,  y  =  x  +  py  -  x  +  6yx  y 

For  u  «  0,  5  =  0,  this  equation  has  the  first  integral 
H(x,y)  =  y2/2  -  x2/2  *  x4/4. 


Some  of  the  level  curves  H(x,y)  =  b  of  this  function  are  shown  in  Fig.  9.3.  For 
b  *  0,  the  curve  is  a  figure  of  eight  and  for  b  >  0,  it  is  a  closed  curve  through 
the  point  ( x,y )  =  (0,  (2b) 1/<2) .  For  b  <  0,  the  set  H(x,y)  =  b  consists  of  two  closed 


Figure  9.3 

curves  surrounding,  respectively,  the  equilibrium  points  (1,0),  (-1,0).  These  curves 

2  4 

pass  through  the  point  (0,c),  0<c<l,b=-c/2*c/4.  The  derivative  H(x,y) 

•  2  2  2 

along  the  solutions  of  (9.2)  is  given  by  H(x,y)  =  py  *  6yx  y  . 

The  first  step  in  the  analysis  is  to  determine  the  curve  p  =  p*(6)  so  that 
(9.2)  has  a  homoclinic  orbit.  This  curve  is  obtained  as  in  the  example  of  Section  8 
and  is  shown  to  be 

U*(6)  -  - (4/5)6y  ♦  0(S2)  as  6-0. 


The  next  step  is  to  analyze  the  periodic  orbits.  There  is  the  possibility  of 
two  types:  an  orbit  which  contains  only  one  equilibrium  point  in  its  interior  or  an 
orbit  which  contains  three  points  in  its  interior.  These  orbits  will  be  close  to  a 
curve  H(x,y)  «  b  for  some  b  <  0  in  the  first  case  and  some  b  >  0  in  the  latter  case. 
These  two  cases  must  be  analyzed  separately. 

For  b  >  0,  let  t.  *  t,(iJ,b,5)  >  0  be  the  first  positive  value  of  t  for 

*  *  i/2 

which  the  solution  through  (0,{2b)  )  crosses  the  x-axis,  say  at  Xj(p,b,6).  Let 

t,  *  T,(u,b,6)  be  the  first  negative  value  of  t  for  which  the  solution  through 
1  1  1/2 

(0,-(2b)  )  crosses  the  x-axis,  say  at  x,(p,b,6).  From  the  symmetry  in  the  equa¬ 
l/2  i 

tion,  it  follows  that  (0,(2b)  )  lies  on  a  periodic  orbit  if  and  only  if 

H(Xj  (v,b,<)  ,0)  *  HfXjfu.b,  S)  ,0) .  Using  the  expression  for  ll(x,y),  it  is  not  diffi¬ 
cult  to  show  that  this  implies 
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u 


V  (b,«)  =  -yP(b)6 

1  ft  2 

x  ydx 


P(b)  = 


f 

Jo 


|Cydx 

J0 


0(62) 


(9.4) 


4  2 

Where  c  =  c(b)  is  the  positive  solution  of  4b  =  c  -  2c  and  H(x,y)  =  b.  For  P(b), 
one  can  now  prove  the  following  basic  result. 


Lemma  9.1.  P(b)  -*•  »  as  b  ♦  ■  and  there  is  a  unique  minimum  of  P(b)  at^  b  =  bj  and 
P"(bj)  >  0. 


Idea  of  the  proof:  If  r(w)  «  (w2-w4/2  ♦  2b)1/l2,  a(b)  » 
4b  -  c4  -  2c2,  then  P(b)  ■  6{b)/a(b).  If  a'  =  da/db.8' 


[  r (w)dw,  8(b)  =  f  w2r(w)dw, 
J0  J0 
=  de/db,  then  one  -shows  that 


3a  =  4ba'  ♦  8' 

158  =  4ba'  +  (4+12b) 8' 


(9.5) 


after  several  computations. 

Now,  suppose  that  P'fbj)  *  0.  Then  a(b1)P"(b1)  *  6" Cbj )  -  Pib^a'^bj)  and 
relations  (9.5)  imply  that 

4b1(4b1*l)(8"(b1)-P(b1)a"(b1)]  =  a'(b1)[P2(b1)*8b1P(b1)-4b1]. 

Thus,  P"(bj)  has  the  same  sign  as  P2 (bj)+8b1P(b1)-4bj . 

On  the  other  hand,  P'(bj)*  0  and  relations  (9.5)  imply  that 

SP2(b1)+8bjP(b1)-4P(b1)-4b1  *  0 

which  implies  P(bj)  <  1  since  bj  >  0.  Using  the  fact  that  SbjPfbj)  -  4bj  * 

4P(bj)  -  5P2(bj),  we  see  that  P"(bj)  has  the  same  sign  as  P(bj)  -  P2 (b ^  which  is 
>  0  since  P(bj)  <  1 . 

It  is  not  difficult  to  prove  that  P(b)  -*  ■»  as  b  +  •  and  P '  (b)  ♦  as  b  +  0. 
This  will  complete  the  proof  of  the  lemma. 

The  graph  of  P(b)  is  illustrated  in  Fig.  9.4. 


Figure  9.4 
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From  (9.4),  the  periodic  orbits  with  b  >  0  are  given  approximately  by  -u(yS)  - 
P(b).  Thus,  approximately,  there  should  be  no  periodic  orbit  encircling  the  three 
equilibrium  points  if  -m(yS)  *  <  min  P(b)  *  P(b.),  two  periodic  orbits  if  P(b  )  < 

-1  1  _i  1 

-u(Yi)  <  P(0)  and  one  periodic  orbit  if  P(0)  <  -u(y6)  .  This  can  be  made  precise 

since  P"(bj)  >  0  to  confirm  the  part  of  the  bifurcation  diagram  in  Fig.  9.1  for  the 
periodic  orbits  which  encircle  three  equilibrium  points. 

The  analysis  of  the  periodic  orbits  encircling  only  one  equilibrium  point  uses 
methods  very  similar  to  the  ones  in  Section  8  and  will  not  be  given. 


These  notes  are  intended  to  relate  the  results  stated  in  previous  sections  to 
existing  literature.  No  claim  is  made  toward  completeness  nor  even  original  sources. 

Section  1.  The  methods  in  this  section  are  very  special  cases  of  a  much  more  gener¬ 
al  global  procedure  for  discussing  the  zeros  of  functions.  This  procedure  often  is 
called  the  alternative  method  and  originated  from  some  fundamental  papers  of  Cesari 
in  the  early  1960's  (for  references  and  an  historical  discussion,  see  Cesari  [6], 
Chow  and  Hale  [7]). 

Section  3.  Theorem  3.1  can  be  found  in  deOliveira  and  Hale  [13]  and  can  also  be 
obtained  from  a  result  in  Golubitsky  and  Schaeffer  [15].  Theorem  3.2  appeared  in 
the  paper  of  Crandall  and  Rabinowitz  [11],  [12].  The  saddle-node  and  cusp  bifurca¬ 
tions  can  be  found  in  Andronov  et  al  [1],  The  complete  discussion  of  the  case 
G(x,0)  *  8xq  ♦  o(|x|q),  8  #  0,  belongs  to  the  general  theory  of  unfolding  of  singu¬ 
larities  (see,  for  example,  Golubitsky  and  Buillemin  [14]).  For  a  full  treatment 
of  the  evolutionary  equations  of  the  form  (3.17),  see  Henry  [17].  For  functional 
differential  equations,  see  Hale  [16]. 

Section  4,  For  a  detailed  discussion  of  the  Hopf  Bifurcation  Theorem  as  well  as 
references,  see  Marsden  and  McCracken  [21]  and  Chow  and  Hale  [7].  Results  and 
references  for  parabolic  equations  may  be  found  in  Kielhbfer  [20]  and  for  functional 
differential  equations  in  [16]. 

Section  5.  An  excellent  discussion  of  dichotomies  is  contained  in  Coppel  [10]. 
Exponential  dichotomies  for  parabolic  equations  are  contained  in  Henry  [17]  and  for 
functional  differential  equations  Pecelli  [24].  Lemmas  5.3  and  5.4  are  due  to 
Palmer  [23].  The  use  of  the  Fredholm  alternative  to  discuss  Example  (5.11)  was 
first  given  by  Chow,  Hale  and  Mallet-Paret  [8].  The  Mel'nikov  function  can  also  be 
used  to  discuss  (5.11)  (see,  e.g.  Holmes  [18]).  Lemma  5.4  and  its  application  to 
obtain  the  function  G(a,u,a)  in  (5.10)  can  be  considered  as  a  generalization  of  the 
Mel’nikov  function  to  n-dimensions. 


Section  6.  The  methods  and  results  in  this  section  are  based  on  Palmer  [23].  For 
other  references  and  approaches  to  the  "shadowing  lemma"  and  the  symbolic  dynamics 
of  Corollary  6.5,  see  Smale  [27],  Conley  [9],  Moser  [22],  Sil’nikov  [26]. 

Section  7.  Theorem  7.1  is  due  to  Andronov  and  Pontrjagin  [2]  and  Peixoto  [25]. 
Theorem  7.2  is  due  to  Andronov  et  al  [1]  and  Sotomayor  [28].  Theorem  7.3  is  due  to 
Sotomayor  [28]. 

Section  8.  The  bifurcation  diagram  for  Eq.  (8.2)  was  considered  by  Howard  and 

Koppell  [19].  Arnol'd  [3]  and  Bogdanov  [4]  considered  an  equivalent  equation 
•  .  2 

x  «  y,  y  *  e-  ♦  e.x  ♦  ax  *  8xy.  Bogdanov  [4]  has  shown  that  every  two  parameter 

^  *  .  .  2  3 

family  of  vector  fields  close  to  x  =  y,  y  =  ax  +  6xy  in  the  C  topology  is  equiva¬ 
lent  to  a  member  of  the  above  two  parameter  family. 

Section  9.  The  results  in  this  section  are  due  to  Takens  [29,30]  and  Carr  [5]. 


Acknowledgement :  The  author  is  indebted  to  Giorgio  Fusco  for  several  stimulating 
discussions  during  the  preparation  of  these  notes. 
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